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1.0  INTRODUCTION 

Considerable  progress  has  been  made  in  the  past  few  years  in  applying  the 
finite  element  method  to  the  analysis  of  shell  structures.  The  initial  work  in  this  area 
was  devoted  to  shells  of  revolution  in  which  closed  rings  or  conical  shell  segments  were 
used  to  model  complete  structures.  Much  of  this  work  is  reviewed  by  Jones  andStrome. 

Attempts  to  develop  a  finite  element  method  for  general  shell  structures  have 
generally  followed  two  different  courses.  In  the  first  approach,  the  shell  is  replaced 
by  an  assemblage  of  flat  plate  elements  that  are  either  triangular  or  quadrilateral  in 
shape.  Each  plate  element  is  connected  in  some  fashion  to  those  surrounding  it,  and 
undergoes  both  bending  and  stretching  deformations.  This  approach  has  been  success¬ 
fully  employed  for  cylindrical  geometry  by  Hrennikoff  and  Tezcan2’  and  for  general  shell 
shapes  by  Zienkicwicz  and  Cheung3’,  Clough  and  Johnson4’,  and  Carr5’.  However,  the 
method  has  the  disadvantage  that  there  is  no  coupling  between  bending  and  stretching 
within  each  element,  and  consequently  a  large  number  of  elements  must  be  used  to 
achieve  satisfactory  accuracy. 

The  second  approach,  which  ultimately  should  yield  better  results,  is  to  de¬ 
velop  curved  shell  elements  that  permit  closer  geometrical  representations  of  a  shell 
structure.  Such  an  approach  has  been  followed  with  good  results  for  the  case  of  cylin¬ 
drical  shells  where  rectangular  elements  are  generally  adequate.  Successful  cylindri¬ 
cal  shell  elements  have  been  developed  by  Bogner,  Fox  and  Schmit6’,  Cantin  and  Clough7, 
and  Olson  and  Lindberg8’. 

On  the  other  hand,  attempts  to  develop  rectangular  or  quadrilateral  elements 
for  general  shells  have  been  only  partially  successful.  For  example,  Connor  and 
Brebbia9’  have  introduced  a  rectangular  element  based  on  shallow  shell  theory,  and 
Gallagher'0’  has  developed  a  quadrilateral  element  for  translational  shells.  However, 
both  elements  use  only  linear  distributions  for  the  tangential  displacements  and  conse¬ 
quently  do  not  incorporate  all  the  required  rigid  body  modes.  Although  this  does  not 
preclude  convergence  to  the  correct  answer,  it  does  mean  that  the  elements  are  far  too 
stiff  and  therefore  are  quite  inefficient. 

Several  curved  shallow  shell  elements  of  arbitrary  triangular  shap^  Mvc  now 
been  developed.  Utku  and  Melosh"’  have  introduced  a  shallow  curved  triangular  element 
based  only  on  linear  displacement  distributions.  This  element  is  consequently  too  stiff. 
More  recently,  Strickland  and  Lodcn12’  have  presented  another  shallow  curved  triangular 
element  based  on  a  cubic  variation  for  the  normal  displacement  and  linear  variations 
for  the  tangential  displacements.  This  element  is  a  distinct  advancement  over  the 
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FOREWORD 


The  Quarterly  Bulletin  is  designed  primarily  for  the 
information  of  Canadian  industry,  Universities,  and  Government 
Departments  and  agencies.  It  provides  a  regular  review  of  the 
interests  and  current  activities  of  two  Divisions  of  the  National 
Research  Council  of  Canada: 

The  Division  of  Mechanical  Engineering 
The  National  Aeronautical  Establishment 

Some  of  thcworkof  the  two  Divisions  comprises  classified 
projects  that  may  not  be  freely  reported  and  contractual  projects  of 
limited  general  interest.  Otherwork,  not  generally  reported  herein, 
includes  calibrations,  routine  analysis  and  the  testingof  proprietary 
products. 


Comments  or  enquiries  relating  to  any  matter  published 
in  this  Bulletin  should  be  addressed  to  the  Editor,  DME/NAE 
Bulletin,  National  Research  Council  of  Canada,  Ottawa,  mentioning 
the  number  of  the  Bulletin. 
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former  one,  but  appears  to  oiler  imle  improvement  in  accuracy  over  the  comparable 
flat  plate  model  of  Clough  and  Johnson41.  Bonnes,  Dhatt,  Giroux,  and  Robichaud ,3!  have 
introduced  a  shallow  curved  triangular  clement  based  on  a  cubic  variation  for  all  three 
displacement  components.  They  actually  present  results  for  two  versions  of  this  ele¬ 
ment,  one  with  27  degrees  of  freedom  and  the  other  with  36.  The  latter  version  incor¬ 
porates  the  normal  derivatives  of  u,  v,  and  w  at  mid-edge  nodes  as  the  extra  degrees 
of  freedom  and,  hence,  allows  the  normal  derivatives  to  vary  quadratically  along  an 
edge,  whereas  the  former  one  has  only  a  linear  variation.  Both  elements  appear  to 
give  somewhat  better  results  than  the  element  of  Strickland  and  Loden'3’.  However, 
the  36-degree-of-freedom  version  has  the  distinct  disadvantage  of  having  mid-edge 
nodes,  and  the  other  model  is  somewhat  limited  by  having  only  a  linear  variation  of 
normal  derivatives  along  each  edge. 

All  the  foregoing  elements  have  disadvantages  common  to  triangular  plate¬ 
bending  elements  as  well,  and  it  is  not  surprising  that  the  development  of  a  good  trian¬ 
gular  curved  shell  element  has  been  impeded  by  the  lack  of  a  satisfactory  conforming 
triangular  plate-bending  element.  Only  recently  has  such  a  plate-bending  element  been 
developed  by  Cowper,  Kosko,  Lindberg  and  Olson, u"'6)  and  also  by  Bell,7'I8)  and  Butlin 
and  Ford'9'.  This  element  uses  as  generalized  co-ordinates  the  transverse  displace¬ 
ment  and  its  first  and  second  derivatives  at  each  vertex,  a  total  of  18  in  all.  The  dis¬ 
placement  function  for  the  element  contains  a  complete  quartic  polynomial  plus  some 
higher  degree  terms  and  allows  a  cubic  variation  of  the  normal  slope  along  each  edge. 

It  is  shown15  ,6>  that  use  of  this  element  leads  to  strain  energy'  convergence  rates  ap¬ 
proaching  n'6,  where  n  is  the  number  of  elements  per  side  of  a  plate.  Along  with  this 
rapid  rate  of  convergence,  it  was  also  found  that  remarkably  accurate  displacement 
and  stress  predictions  were  obtained,  even  with  coarse  grids  of  elements. 

The  following  article  is  concerned  with  extending  this  highly  successful  ap¬ 
proach  to  the  analysis  of  arbitrarily  shaped  shell  structures.  As  a  first  step  in  this 
direction,  the  application  to  shallow  shells  is  investigated,  and  a  shallow  shell  element 
of  general  triangular  shape  is  formulated  and  tested.  This  element  proves  to  be  excep¬ 
tionally  efficient,  yielding  displacement  and  stress  predictions  of  comparable  accuracy 
to  that  of  the  plate-bending  element  discussed  above.  Finally,  the  methodology  for 
applying  this  shallow  shell  element  to  arbitrary  deep  shells  is  developed  and  tested  with 
numerical  applications. 


2.0  A  REFINED  SHALLOW  SHELL  FINITE  ELEMENT 

In  this  section,  a  brief  derivation  of  the  refined  shallow  shell  element  and 
numerical  results  from  one  application  are  given.  A  more  complete  development  and 
additional  numerical  results  arc  available  in  Reference  20. 

2.1  Theoretical  Formulation 

2.1.1  Strain  Energy 

The  geometry  for  an  arbitrary  triangular  shallow  shell  element  is  shown  in 
Figure  1.  The  shell  shape  is  defined  by  the  height  £(4,tj)  above  the  base  plane,  in  which 
£,T7  are  taken  as  local  co-ordinates  and  x,y  as  global  co-ordinates.  The  dimensions 
a,  b,  c  of  the  base  triangle  1  ’ 2* 3’ ,  and  the  rotation  angle  0,  are  easily  derived  in  terms 
of  the  global  co-ordinates  of  the  vertices15  ,6) 
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l'ollowing  the  shallow  shell  theory  of  Novozhilov,2"  the  membrane  strains  in 
the  shell  arc  given  by 


u 

ut  - 

£u  w 

in 

V*  - 

£i.i,w 

(1) 

■fl 

vf  “  2£{„  w 

where  u  and  v  are  the  tangential  displacements  measured  parallel  to  the  £  and  rj  di¬ 
rections,  respectively,  and  w  is  the  normal  displacement  (Fig.  1).  The  subscripts  on 
u,  v,  w  and  t  denote  differentiation,  i.e.  u{  -  8u/8|  etc.  The  bending  strains  are 
given  by 
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Combining  the  contributions  from  the  membrane  and  bending  strains  yields  the  strain 
energy  density 
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For  shallow  shells,  the  area  of  the  shell  surface  is  approximately  equal  to  its  projected 
area,  and  hence  equation  (3)  may  be  integrated  over  the  base  plane.  Then  combining 
equations  (1)  to  (3)  yields  the  final  strain  energy  expression 
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where  the  integration  is  over  the  base  triangle  1'2'3'  in  Figure  1. 

In  the  present  work,  the  function  is  assumed  to  be  of  quadratic  form 
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rhis  implies  that  the  shell  element  has  constant  curvatures,  and  this  is  consistent  with 
the  approximations  oi  shallow  shell  theory.  It  may  be  noted  from  equation  (4)  that  the 
sneli  curvatures  are  the  only  shape  quantities  required  in  calculating  strain  energy, 
therefore,  they  may  be  specified  for  an  element  directly  (which  is  the  approach  followed 
n  this  bection),  or  they  may  be  calculated  from  equation  (5)  once  the  constants  therein 
have  been  determined.  One  convenient  way  to  do  this  is  to  specify  the  height  of  the 
shell  at  points  1',  2',  and  3'.  and  also  at  the  mid-edge  points  of  the  element  of  Figure  1, 
ihis  is  the  approach  followed  in  the  deep  shell  applications  of  the  next  Section. 

2.1.2  Displacements 

Special  displacement  functions  are  required  for  the  present  element  in  order 
to  ensure  conformity  and  high  accuracy.  The  displacement  function  developed  for  the 
conforming  plate-bending  element  (lief,  15,  13)  is  used  here  for  the  normal  displace¬ 
ment  w.  In  Reference  20,  it  is  shown  that  the  use  of  complete  cubic  polynomials  for 
each  of  the  tangential  displacements  u  and  v  leads  to  an  asymptotic  strain  energy  con¬ 
vergence  rate  of  n'A,  where  n  is  the  number  of  elements  per  side  of  a  shallow  shell. 

Hcnv'i,  the  starting  point  for  the  present  shallow  shell  element  is  to  assume 
u,  v,  ana  w  in  t,  c  form  (Fig.  1) 

u  =  m  +  a2(;  +  i>v  +  a4$7  *-  a34r?  +  a6y2  +  a743  u  a84J»l  *  a&l2  +  a,0n3  (6a) 

v  =  an  i  n»{  +  an.*  +  au4J  ••  al34rj  +  a,6y3  +  a, ,4 3  +  nw4*»J  +  +  a?cfl3 

(6b) 

w  =  a2t  +  h  W1  +  a  it,1  +  a234n  +  a j6tj2  +  a7743  +  a2a4‘TJ 

+  a294hJ  ajoi)3  +  a3i44  +  a3243r)  +  a334V  +  n344»?3 
+  a35v/4  +■  a3643  +  a374V  +  a3S4V  +  a  394';4  +  a40rj 5  (6c) 


Note  that  the  expression  for  w  (eq.  (6c))  docs  not  contain  the  term  4  V,  and  therefore 
automatically  satisfies  the  requirement  that  the  normal  slope  be  only  cubic  along  the 
clement  edge  1-2.  The  two  conditions  that  ensure  only  cubic  variations  of  normal 
slopes  along  edges  2-3  and  3-1  arc  the  same  as  for  the  flat  plate  casc,5,,6’30>.  These 
constraints  are  just  sufficient  to  reduce  the  twenty  independent  parameters  of  equa¬ 
tion  (6c)  to  only  eighteen.  The  generalized  displacements  for  w  are  then  taken  to  be 
w,  and  its  first  and  second  derivatives  at  each  corner  of  the  clement,  a  total  of  eighteen, 
which  is  consistent  with  the  eighteen  free  parameters  available.  The  generalized  dis¬ 
placements  for  u  and  v  arc  taken  to  be  u  and  v  and  their  first  derivatives  at  each 
corner  plus  u  and  v  at  the  element’s  centroid.  This  gives  a  total  of  twenty,  which  is 
consistent  with  the  twenty  free  parameters  of  equations  (6a)  and  (Gb),  All  the  general¬ 
ized  displacements  arc  assembled  into  a  38-column  vector  Jw,j,  firot  in  the  local  co¬ 
ordinate  system  4,7? 
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where  u{  »u/l>4,  etc.  The  subscripts  1,  2,  3,  c  denote  the  corners  1,  2,  and  3  and 
Hie  centroid  oi  the  element,  respectively.  The  coefficients  a.  of  equations  (6)  are  as¬ 
sembled  into  a  40-column  vector  J  A  {  where  ' 


f A i  -  (a],  3j,  .  .  ,  ,  ,  u^q)  (8) 


and  combining  equations  (6)  to  (8)  plus  the  two  constraint  equations  yields  the  matrix 
relation 


(W, 

/  0 
\  0 


(9) 


where  the  40  y  40  transformation  matrix  [t]  is  given  in  Table  I, 

The  determinant  of  [t]  has  the  value  -64  c34  (a  +  b)31  (a2  +  e2)  (b3  +  c2)/729, 
which  is  nonzero  for  all  practical  problems.  Hence,  equation  (9)  may  be  inverted  to 
gi  vc 


JA|  -  [t-'J  j  0  |  =  [t,]  }\v, !  (10) 

where  the  40  x  38  matrLx  [ T , ]  consists  of  the  first  38  columns  of  [l’"']. 

2.1.3  Stiffness  MatrLx 


The  stiffness  matrLx  for  the  element  is  obtained  from  a  calculation  of  strain 
energy.  The  displacement  functions  of  equations  (6)  arc  substituted  in  equation  (4)  and 
the  integration  carried  out  to  yield  the  quadratic  strain  energy  form 


u  a  2U-U  :a!  [k]!At 

(11) 

Tile  entries  of  the  stiffness  matrLx  [k]  may  be  determined  in  closed  form  just  as  they 
were  for  the  plate-bending  element13  t*’.  This  may  be  carried  out  most  easily  by  re¬ 
writing  equations  (G)  as 
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Then  substituting  equations  1 12)  into  equation  (4;  and  incorporating  the  sjmmi  try  re¬ 
quirement  yields 
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-  (£„„  +  v  cu)  [q,F  (p,  +  r  ,  q,  +  s,  -  1)  +  q,F  (P|  +  r, ,  q;  +  s,  -  1)] 

(13) 

-  (I-*')  £{„  [n,F  (rn,  *  r  ,  n,+  sj  -  1)  +  n  F  (mj  *  rlt  u(  +  s,  -  1) 

+  P,F  (P,  +  r  -  1,  q,  +  sj)  +  p(F  (Pj  +  r,  -  1,  q(  +  s,)] 

+  fo/  +  £m>*  +  2v  £((£vv+  2(1-^)  Un2)  F  <r.  +  r,>  s.  +  s,> 

*  (t2/ 12)  ]  r, r(  (r,  -  1)  (rj  -  1)  F  (r,  +  r,  -  4,  s,  +  S() 

+  s.s,  (s,  -  1)  (SJ  -  1)  F  (r,  +  rJt  s,  +  sj  -  4) 

+  [2(l-v)  r, Tj s, s,  +  v  r.s,  (r,  -  1)  (sj  -  1) 

+  v  r J s ,  (r  -  1)  (s,  -  1)]  F  (r,.  +  r,  -  2,  s,  +  s)  -  2)| 


where 


F  (m,  n)  =  cn  +  1  [am+ 1 


mi  ni 

(m  +  n  +  2)i 


(14) 


Note  that  i  and  j  run  from  1  to  40,  and  theretore  m,  and  n,  are  defined  to  be  zero  for 
i  >  10,  p,  and  q,  are  zero  for  i  <  11  and  i  >  20,  and  r,  and  s,  are  zero  for  i  <  21.  All 
the  computations  involved  in  evaluating  the  k  from  equation  (13)  are  carried  out  within 


the  computt  r  unit  the  values?  of  a  b,  e,  m  ,  n  ,  p  ,  qr,  r  ,  and  s.  are  furnished,. 


The  transformation  from  local  to  global  co-ordinates  is 


(15; 


where 


)\V2;  -  (U„  u,...  Uy„  V,,  V, v  „  w.„  wyl,  w„,,  wly!, 


w„  U2, ....  u3,,..,  uc,  vt) 


(16) 


is  the  generalized  displacement  vector  relative  to  global  co-ordinates  (u,  v  in  the  direc¬ 
tions  of  x,  y),  and  [r]  is  the  rotation  matrix  given  in  Table  II.  Then  combining  equa¬ 
tions  (10),  (11),  ana  (15)  yields  the  strain  energy  in  terms  of  [W2(  as 


U 


Et 

2(1-1/ 2) 


w,; 


where 


K)  *  [R]'[T,l'lk]lTl]lR] 


(17; 


(18) 


is  the  38  x  38  stiffness  matrix  relative  to  global  co-ordinates. 

2.:1.4  Consistent  Load  Vector 

The  consistent  load  vector  is  obtained  by  calculating  the  virtual  work  done  by 
the  applied  loads  Q^.rj),  Q v(4,?]),  andQw(4,r;)  in  the  u,  v,  and  w  directions,  respec¬ 
tively.  The  transformed  load  vector  becomes 

!p{  =  Mb,]  Iq(  as) 


where  the  entries  in  the  column  vector  j  Q  {  are 


II  $"V'  d$  dTL  1  <  n 

Q,  =)  II  QV4P|  7?q'  *1,  io  <  i  <  21 


\  If  QJ''  0s'd4  drj,  i  >  20 


(20) 
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2.1.5  Condensation  of  Stiffness  Matrix 

Before  proceeding  to  use  the  element  just  derived,  it  is  convenient  to  con¬ 
dense  out  the  centroidal  displacements  u£  and  vc.  It  may  be  noted  that  these  displace¬ 
ments,  since  they  lie  inside  the  element,  will  be  unaffected  when  the  elements  are 
joined  together  to  represent  a  structure.  Hence,  they  may  be  eliminated  before  the 
elements  are  joined  without  affecting  the  final  results.  The  reduction  is  carried  out  by 
minimizing  the  potential  energy  in  one  element.  The  matrix  equation  of  equilibrium 
for  an  element  is  written  in  the  partitioned  form 


where  jwc(  is  the  two-component  vector  (uc,  vc),  and  jw{  contains  the  first  36  com¬ 
ponents  of  equation  (16).  Equation  (21)  is  separated  into  two  equations  and  )  Wc  |  is  elim¬ 
inated  to  yield 

[k])w(  =  )p[  (22) 

where 

[k]=  [K0]  -  [koc][kc  ,][koc]T  (23) 

is  the  reduced  36  x  36  stiffness  matrix,  and 

}P(  =  )Po(  -  [Kot][Kc'’])Pc*  <24) 


is  the  reduced  36  consistent  load  vector  for  the  shell  element. 
2. 2  Numerical  Application  -  Cylindrical  Shell  Roof  Problem 


The  geometry  for  this  problem  is  shown  in  Figure  2,  The  shell  is  loaded  by 
its  own  dead  weight  and  is  supported  by  diaphragms  at  the  ends  but  is  free  along  the 
sides.  Using  symmetry,  only  one-quarter  of  the  shell  was  analyzed,  and  various  uni¬ 
form  gridworks  of  elements  were  used  with  orientations  as  shown.  The  numerical  data 
given  in  the  Figure  were  used  in  the  calculations,  so  that  the  results  could  be  compared 
directly  with  those  of  other  authors. 

Some  of  the  most  pertinent  numerical  results  are  given  in  Table  III  and  con¬ 
vergence  plots  for  these  quantities  are  given  in  Ffgure  3 .  It  is  seen  that  the  approxi¬ 
mations  otU-’ned  with  even  the  coarsest  gridworks  are  exceedingly  accurate,  and  the 
convergence  is  very  rapid.  In  particular,  the  slopes  of  the  error  curves  for  displace¬ 
ments  and  strain  energy  are  all  about  -5. 
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The  quantity  that  has  been  used  by  many  authors  for  comparison  purposes  is 
the  vertical  deflection  of  point  B  (Fig.  2)  .  Such  a  comparison  is  given  in  Figure  4, 
where  the  abscissa  is  the  total  number  of  degrees  of  freedom  (before  boundary  condi¬ 
tions)  required  with  the  various  finite  element  representations.  It  is  clear  that  this 
comparison  favours  the  present  results  in  that  it  shows  that  the  present  element  yields 
the  best  accuracy  for  the  minimum  number  of  degrees  of  freedom.  It  is  interesting  to 
note  that  the  results  of  References  5  and  13  have  not  actually  converged  to  the  exact 
solution,  even  though  more  than  1000  degrees  of  freedom  were  used. 


3.0  ANALYSIS  OF  ARBITRARY  DEEP  SHELLS 

In  this  Section,  a  brief  description  is  given  of  the  development  necessary  to 
apply  the  foregoing  shallow  shell  element  to  arbitrary  deep  shells.  This  is  followed  by 
some  results  from  two  numerical  examples.  All  this  work  is  preliminary  in  nature, 
but  full  details  will  be  presented  in  an  NRC  Aeronautical  Report  when  the  investigations 
are  completed. 

3.1  Theoi-etical  Development 

3. 1. 1  Co-Ordinate  Systems 

The  relevant  geometry  is  shown  in  Figure  5,  where  X,  Y,  Z  are  cartesian 
global  co-ordinates,  4,rj,f  are  cartesian  local  co-ordinates,  and  a, /3  are  curvilinear 
shell  co-ordinates  lying  in  the  mid-surface  of  the  shell.  At  this  point,  the  latter  need 
not  be  orthogonal.  The  4,tj  co-ordinates  are  to  be  used  as  the  base  plane  co-ordinates 
in  the  shallow  shell  formulation  and  are  defined  to  go  through  the  element  corner  nodes 
1,  2,  3  as  shown. 

Now  it  is  explicitly  assumed  that  a  point  a,p  on  the  shell  may  be  determined 
in  terms  of  the  global  co-ordinates  by 


X  =  X  (a,/?) 

Y  =  Y  (<y,/3)  (25) 

Z  =  Z  {a, 13) 


This  leads  to  a  natural  computation  process  in  which  the  element  corner  nodes  are 
located  first  by  specifying  their  shell  co-ordinates  a, (3,  and  then  their  global  co¬ 
ordinates  X,  Y,  Z  are  calculated  from  equations  (25).  The  element  geometry  may  then 
be  determined  from  these  latter  co-ordinates. 

Some  results  from  analytical  geometry  follow.  The  global  co-ordinates  of 
the  £,tj.£  system  origin  arc 

x0  =  (1-p)  X,  +  pX2 
Y0  =  (1-p)  Y,  +  pY j 
Z0  =  (1-p)  Z,  +  pZ2 


(26) 
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where  the  parameter  p  is 


-(X.?  r  M  (X?  -  X,)  +  (Y3  -  Y|)  (Y?  -  Y,)  +  (Z3-  Zj)  (Z?  -  Z.) 
(Xj-  X,)2  *•  (Y:  -  Y))z  +  (Z2-  Z,)» 


(27) 


The  subscripts  denote  the  element  corner  nodes  1,  2,  3  of  Figure  5.  The  element  dimen¬ 
sions  a,  b,  c  are  then  given  by 

a  =  (1-p)  /(Xj  -  X,)?  +  (Yj  -  Y,)J  >  (Z2-  z7 

b  =  P  J  (X,-  X,)2  +  ( Yj  -  Y,)2  +  (Z2  -  Z,7  (2«) 

c  =  ^(X3-Xo)2  +  (Y3-Y0)j  +  (Z3  -  Z0)? 


3.1.2  Relations  Between  Co-Ordinate  Systems 

The  following  relations  between  the  global  and  local  cartesian  co-ordinates 
are  obtained 


4  =  [(X2  -  Xt)  (X  -  XJ  +  <Y,  -  Y,«  (Y  -  Y0)  +  (Z2  -  Z,)  (Z  -  Z0)]/(a  +  b) 


V  -  [(X3  -  x0)  (X  -  XJ  +  (Y3  -  Y0)  (Y  -  Y0)  f  (Z3  -  Z0)  (Z  -  Z0)]/c  (20) 

r.  ■---  b,  (X  -  X0)  +  b,  (Y  -  Y0)  +  b3  (Z  -  Z0) 


where 


b,  =  [(Y,  -  Y,)  (Z3  -  Z0)  -  (Z,  -  Z|)  <Y3  -  Y0)]/(a  +  b)  c 

b2  =  [(Z2  -  Z,)  (X3  -  X0)  -  (X2-  X,)  (Z3  -  Z0)]/(a  +  b)  c  (30) 

b3  -  [(Xj  -  X,)  (Y3  -  Y0)  -  (Yj  -  Y,)  (X3  -  Xo)]/(a  +  b)  c 

Combining  equations  (25)  and  (29)  yields  the  base  plane  co-ordinates  as  functions  of  the 

shell  co-ordinates,  written  symbolically  as 


«  =  «  («,P) 

V  ~  V  (a,p) 


(31) 
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Now  ali  the  derivatives  required  in  transforming  the  generalized  displace  men* c  from 
local  co-ordinates  to  shell  co-ordinates  may  be  derived  from  the  above  equations.  Un¬ 
fortunately,  these  derivations  soon  become  unwieldly,  and  hence  a  simplifying  assump¬ 
tion  is  introduced  here  to  reduce  the  complexity.  As  a  first  approximation,  it  is  as¬ 
sumed  that  equations  (31)  may  be  expanded  to  first  order  about  each  element  node,  say 


$  ®  d,  -  e,,a  -  e)it? 
rj  =  dj  *  e,,a  -  e-^ 


(32) 


where 


X2  -  Xi  ax  |  Y2  -  Yi  3Y 

Z2  -  Zi  az 

a  ^  b  da  a  b  da 

a  +  b  da 

Xi  -  Xi  3X  (  Y2  -  Yi  9Y 

72  -  Zi  az 

a  -*■  b  3/3  a  *-  b  9 /3 

a  L  b  9/3 

X3  -  Xo  ax  1  Y3  -  Y0  9Y 

z3  -  Zo  a  z 

c  da  c  da 

c  act 

X3  -  Xo  ax  Yj  -  Yo  9Y 

Z3  -  Zo  9  Z 

c  dp  0  9/3 

c  dp 

(33) 


are  to  be  evaluated  at  an  element  node.  Equations  (32)  will  be  valid  approximations 
within  seme  neighbourhood  of  this  node  provided  there  are  nc  singularities  there.  As 
the  finite  element  size  is  reduced,  these  neighbourhoods  around  each  node  will  eventu¬ 
ally  overlap,  and  hence  provide  a  uniformly  valid  approximation  over  the  entire  element. 

In  the  following  derivations,  it  will  be  convenient  to  have  the  inverse  of 
equations  (32),  say 

a  =  f,  +  g,jT) 

P  =  *f2  +  *  g22r> 

where 

gn  =  e22/e 

g  j  2  =  -eu/e 

62 1  =  ~e2i/ e 

622  ~  e 1 1/ e 


(34) 


(35) 


and 


e 


elle22  "  el2e2l 
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The  shell  element  curvatures  are  determined  locally  by  specifying  the  shell 
elevations  above  the  £,tj  base  plane  at  the  mid-side  nodes  4,5,6  (Fig..  5)  and  using  them 
in  the  quadratic  function,  equation  (5).  Note  that  the  elevations  are  zero  at  the  corner 
nodes  1,  2,  3..  The  locations  of  the  mid-side  nodes  are  determined  first  in  the  shell  co¬ 
ordinates  simply  as  averages  of  the  corner  node  co-ordinates,  i.e. 


o4  -  (a,  +  a2)/  2, 
a5  =  (n2  -  a3)/2, 

a6  =  +  aa)/2. 


=  (0,  *P2)/  2 
+/33)/2 
h  =  +02)/2 


(36) 


The  global  co-ordinates  of  these  nodes  are  then  easily  obtained  from  equations  (25)  and 
are  substituted  into  equations  (29)  to  yield  the  three  elevations  £  |  and  the  corresponding 
base  plane  co-ordinates  |  ,  rji  for  i  -  4,5,6.  Note  that  these  points  do  not  necessarily 
coincide  with  the  mid-side  points  of  the  base  plane  triangle.  These  results,  together 
with  the  zero  elevations  at  the  three  corner  nodes,  are  then  sufficient  to  complete  the 
quadratic  fit  to  the  shell  shape  and  determine  the  c,  coefficients  of  equation  (5).  Fi¬ 
nally,  the  three  shell  curvatures  required  in  the  stram  energy  calculation  are  obtained 
simply  by  differentiating  equation  (5). 

3.1.3  Transformation  Matrices 

The  shallow  shell  formulation  of  Section  2.  0  begins  with  generalized  dis¬ 
placements  written  relative  to  the  base  plane  co-ordinates,  w,  w£ ,  w,,  etc.  To  use 
these  results  in  the  present  formulation,  it  is  merely  necessary  to  transform  these  base 
plane  generalized  displacements  to  shell  co-ordinate  ones.  The  displacement  normal 
to  the  shell  w  is  a  scalar  and  hence,  using  equations  (34),  its  derivatives  transform 
simply  as 


w  * 

w 

w£ 

Wa 

w 

•  M 

i" 

w  aa 

\  ww 

where  [r6]  is  given  in  Table  IV.  The  linear  approximation  (34)  implies  that  the  second 
derivatives  of  a,/3  with  respect  to  £,rj  are  neglected  in  the  transformation  (37).  The 
tangential  displacements  u,v  are  defined  to  be  parallel  to  the  co-ordinate  axes  and 
hence  must  be  transformed  as  vectors.  Again  using  results  from  analytical  geom¬ 
etry  and  the  element's  shallowness,  together  with  equations  (34)  for  derivatives, 


i 

i 
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leads  to  the  approximation 


\ 

i 

U  1 

u 

i 

i 

U  a 

\  Uv 

i 

■  N  • 

IT1, 

r 

1: 

i 

i 

u  1 

v  r. 

(C8> 


where  u  and  v  are  in  the  directions  of  o  and  and  [ H 5  ]  is  giver.  in  Table  IV,.  Now 
combining  all  these  results  and  dropping  the  "tilda"  notation  for  simplicity,  the  com¬ 
plete  transformation  for  one  element  is 


W 


3* 


(39) 


where 


u0„  ’J/n,  v,. 

v„i, 

V<j!>  Wl>  Wctl>  W0l>  Woal> 

1  » 

U2, 

......  u3,  .....  Ut,  Vt ) 

(40) 

is  the  generalized  displacement  vector  in  shell  co-ordinates  (u,  v  in  the  directions  of 
Q,  (3)  and  [rt4]is  the  new  "rotation"  matrix  given  in  Table  IV „  The  strain  energy  in  the 
shell  element"  then  becomes 


u  “  IfTT*)  <41> 

where 

M  -■  [T,]  !.!<][T,][H4]  (42) 

is  he  38  x  38  stiffness  matrix  relative  to  shell  co-ordinates.  Finally,  the  consistent 
load  vector  and  stiffness  matrix  condensation  procedures  (S“e..  2.1.4  and  2.1.5)  follow 
through  just  as  before  with  [l\]  replaced  by  [U4],  etc. 


The  usual  procedure  tor  assembling  individual  element  stiffness  and  load 
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matrices  into  the  master  matrices  for  a  given  finite  element  gridwork  is  followed  here 
as  well.  This  means  that  all  corresponding  generalized  displacements  are  equated  at 
a  node  (except  when  a  discontinuous  variable  is  expected,  e.g.  normal  curvatures 
across  a  step  in  shell  thickness)  and  sums  of  corresponding  generalized  forces  are  set 
equal  to  the  appropriate  load  components  as  well.  The  question  of  what  continuity  the 
former  procedure  provides  here  is  as  yet  unanswered  and  requires  further  investigation. 

3.  2  Numerical  Applications 

Preliminary  results  have  been  obtained  from  applying  the  foregoing  develop¬ 
ments  to  two  numerical  examples.  These  initial  results  are  very  encouraging  and  are 
presented  in  the  following. 

3.2,1  Pinched  Cylindrical  Shell 

The  problem  geometry  illustrated  in  Figure  6  is  that  of  a  thin  circular  cylin¬ 
drical  shell  of  radius  R,  thickness  t  and  length  L.  The  cylinder's  ends  are  freely 
supported  (sometimes  called  diaphragm  ends) ,  and  it  is  loaded  by  two  centrally  located 
and  diametrically  opposed  concentrated  forces  P.  Using  the  double  symmetry  available, 
only  one-eighth  of  the  cylinder  was  analyzed. 

The  obvious  choice  for  shell  co-ordinates  a,/3  are  the  axial  and  circumfer¬ 
ential  ones  indicated  in  the  Figure.  Equal  subdivisions  of  a  from  0  to  L/2  and  0  from 
0  to  ir/2  were  used  to  form  2  x  2,  3  x  3,  4  x  4,  and  5x5  uniform  grids  of  elements; 
the  3x3  grid  is  shown  in  Figure  6.  A  non-uniform  5x5  grid  with  element  width 
ratios  of  1:  2:  3:  4:  5  going  away  from  the  point  load  (also  shown  in  Fig.  6)  was  used  as 
well.  Numerical  calculations  were  carried  out  for  a  shell  with  R/t  =  100,  L/R  =  2, 
and  v  =  0.3.  Some  typical  numerical  results  are  presented  in  Table  V  and  distribu¬ 
tions  are  plotted  in  Figures  7a  to  d  and  8.  The  exact  results  are  from  a  double  Fourier 
series  solution  of  Flttgge's  equation's”’,  using  80  terms  in  each  direction. 

Ignoring  the  non-uniform  grid  results  of  Table  V  for  the  moment,  it  appears 
that  the  other  finite  element  predictions  all  converge  monotonically  towards  the  exact 
solutions.  On  the  other  hand,  this  convergence  is  quite  slow.  For  instance,  a  plot  of 
the  error  inwc,  the  displacement  under  the  point  load,  versus  n  (for  an  n  x  n  grid) 
reveals  a  slope  of  about  -2. 5.  Now  since  wc  is  linearly  related  to  strain  energy,  this 
means  that  its  convergence  rate  is  also  -2.5.  Although  this  convergence  rate  appears 
low,  it  is  a  direct  consequence  of  the  severity  of  the  pointy  load.  Rates  as  low  as  -2 
have  been  observed  in  some  point-loaded  plate  problems 

The  non-uniform  5x5  grid  yields  an  appreciable  improvement  in  accuracy 
in  the  vicinity  of  the  point  load,  i.  e.  at  C.  Even  the  bending  moments  M„  and  M„  are 
increased  in  magnitude  there.  However,  these  improvements  near  C  are  accompanied 
by  some  degradation  of  accuracy  away  from  the  load,  i.e.  at  B. 

The  plotted  distributions  of  Figures  7a  to  d  provide  a  clearer  picture  of  the 
convergence  trends  in  the  calculated  results.  Firstly,  the  displacements  of  Figure  7a 
exhibit  a  simple  improvement  with  grid  refinement.  Next,  the  membrane  stress  and 
bending  and  twisting  moment  distributions  of  Figures  7b  to  d  also  exhibit  noticeable 
improvements  with  grid  refinements.  However,  there  are  a  number  of  noteworthy  fea¬ 
tures  The  kinks  or  slope  discontinuities  at  element  nodes  first  observed  in  Refer¬ 
ence  20  are  again  in  evidence  in  some  of  these  plots.  This  shows  again  the  importance 
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of  calculating  distributions  of  membrane  stresses  and  moments  rather  than  relying  on 
only  the  nodal  values.  It  appears  from  these  figures  that  N„  near  the  point  load  is  the 
most  difficult  quantity  for  the  finite  elements  to  predict.  In  particular,  the  5x5  grid 
prediction  of  along  DC  is  still  far  from  the  exact  result. 

Most  of  the  above  difficulties  are  solved  by  the  non-uniform  5x5  grid. 

Typical  membrane  stress  distributions  obtained  from  this  grid  arc  shown  in  Figure  8. 

The  improvement  in  N0  along  DC  (which  was  the  poorest  prediction  from  the  uniform 
grids)  is  clearly  outstanding,  the  major  error  remaining  being  that  of  the  small  up-turn 
at  C.  Similar  improvements  in  all  the  other  distributions  were  also  obtained  from  the 
non-uniform  5x5  grid. 

3.2.2  Pinched  Spherical  Shell 

The  second  numerical  application  is  the  pinched  spherical  shell  shown  in 
Figure  8.  The  concentrated  loading  at  each  pole  is  transmitted  to  the  shell  through  an 
infinitely  rigid  boss  with  a  half  angle  of  (30  as  shown.  The  shell  is  assumed  rigidly  at¬ 
tached  or  welded  to  each  boss,  and  hence  the  boundary  condition  there  becomes  one  of 
constraining  v,  w,  and  0w/0 (i  to  allow  only  vertical  motion  of  the  boss. 

The  meridional  angle  a  and  eolatitude  0  are  the  logical  choice  for  shell  co¬ 
ordinates  for  this  problem.  Using  the  symmetry  about  the  equator  and  the  poles,  only 
an  angular  slice  of  the  shell  Is  analyzed  (2x5  grid  shown) .  Although  only  one  element 
is  needed  in  the  meridional  direction,  two  elements  are  used  here  to  provide  a  symmetric 
grid  and  hence  symmetric  numerical  results.  The  meridional  width  of  the  gridwork  a0 
used  in  the  calculations  was  7r/(2n+2),  where  n  is  the  number  of  colatitude  elements. 
Calculations  were  carried  out  for  a  shell  with  R/t  -  50,  v  -  0.3,  and  (Je  -  lOdegrees, 
with  various  values  of  n.  An  exact  solution  for  this  problem  was  obtained  for  compari¬ 
son  purposes  by  numerical  integration  of  Timoshenko's  equations 7 a>. 

The  calculated  results  arc  given  in  Table  VI  for  various  displacements,  mem¬ 
brane  stress  resultants,  and  bending  moments  for  values  of  n  of  3,  5,  7,  and  9,  These 
results  must  be  considered  as  only  preliminary,  since  there  arc  several  factors  yet  to 
be  investigated,  However,  they  arc  very  encouraging,  It  is  seen  that  all  these  results 
appear  to  be  converging  monotonically  towards  the  exact  values  as  the  number  of  ele¬ 
ments  Increases,  Further,  when  the  relative  errors  in  the  displacements  at  the  boss 
wg  and  vg  are  plotted  versus  n,  they  exhibit  convergence  proportional  to  about  n*,s. 

This  implies  that  the  vertical  displacement  of  the  boss,  and  therefore  the  strain  energy 
in  the  shell,  converge  as  n'1  3  as  well. 

Distributions  of  the  displacements  v  and  w  as  obtained  from  the  finite  ele¬ 
ments  are  plotted  in  Figure  10.  The  excellent  accuracy  of  these  predictions  is  evident 
from  the  Figure.  At  this  time,  further  work  is  underway  to  develop  a  method  for  calcu¬ 
lating  the  analogous  membrane  stress  resultant  and  bending  moment  distributions.  It 
should  be  noted  that  the  nodal  predictions  of  those  quantities  appear  to  be  quite  accurate. 


4.0  CONCLUDING  REMARKS 

The  finite  clement  analysis  of  shells,  as  outlined  herein,  is  undergoing  fur¬ 
ther  development.  Several  interesting  and  fundamental  questions  that  arc  under  critical 
investigation  are  as  follows.  What  displacement  continuity  does  the  transformation  to 
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deep  shell  co-ordinates  provide?  What  are  the  effects  of  different  displacement  bound¬ 
ary  conditions?  What  gains  may  be  achieved  through  the  use  of  a  more  exact  transfor¬ 
mation  between  the  local  and  shell  co-ordinates?  What  gains  may  be  obtained  from  a 
more  exact  deep  shell  theory  rather  than  shallow  shell  theory? 

The  results  of  this  investigation  indicate  that  accurate  solutions  for  a  large 
variety  of  shell  problems  can  be  efficiently  obtained.  For  shallow  sheds,  the  present 
approach  using  curved  elements  has  proven  to  be  far  superior  to  the  alternate  flat  cle¬ 
ment  approach.  Though  no  comparisons  have  as  yet  been  made,  it  is  anticipated  that 
the  same  conclusions  will  also  hold  for  deep  shell  problems. 
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TABLE  HI 


NUMERICAL  RESULTS  FOR  CYLINDRICAL 
SHELL  ROOF  PROBLEM 
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TABLE  V 

NUMERICAL  RESULTS  FOR  PINCHED  CYLINDRICAL  SHELL 
R/t  =  100,  L/R  =  2,  v  =  0.3 
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TABLE  VI 


NUMERICAL  RESULTS  FOR  PINCHED  SPHERICAL  SHELL 
R/t  =50,  v  =  0.3 

(SUBSCRIPTS,  B  =  BOSS,  E  =  EQUATOR) 
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FIG. 2:  CYLINDRICAL  SHELL  ROOF  CONFIGURATION 

(3x3  GRID  SHOWN) 
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FIG  3-  RELATIVE  ERROR  OF  FINITE  ELEMENT  SOLUTIONS 
FOR  CYLINDRICAL  SHELL  ROOF  PROBLEM 


FIG, 4:  CYLINDRICAL  SHELL  ROOF  PROBLEM  COMPARISONS 
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DISPLACEMENT  DISTRIBUTES  FOR  PINCHED 
CYLINDRICAL  SHELL 
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FIG.  9  -  PINCHED  SPHERICAL  SHELL  CONFIGURATION 
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APPENDIX  A 
NOMENCLATURE 


Symbol 


Definition 


“i 

|A 

a,  b,  c 
E 

[k].  [k] 

L 

^  i »  ^ i >  Pii  Qii  ^"ii  s j 

M<j,  Ma^j) 

Nor i  N^,  N0^j 

P 

R 

[K] 

t 

[t] 

U,  V 

U 

w 

jw| 

x,  y,  i  ) 
x,  y,  z{ 

a,  p 
0 
v 

4.  V,  t 
subscripts 


Displacement  polynomial  coefficients 
Column  vector  of  a. 

Finite  element  dimensions 
Young's  modulus 
Stiffness  matrices 
Shell  length 

Exponents  of  £  and  t)  in  displacement  polynomials 
Customary  moment  and  membrane  stress  resultants 

Magnitude  of  point  load  on  shells 
Shell  radius  of  curvature 
Rotation  or  transformation  matrix 
Shell  thickness 
Transformation  matrix 
Tangential  displacements  of  shell 
Strain  energy 

Normal  displacement  of  shell 

Column  vector  of  generalized  displacements 

Global  cartesian  co-ordinates 

Shell  co-ordinates 

Element  rotation  angle 

Poisson's  ratio 

Local  cartesian  co-ordinates 

Denote  differentiation,  c.g.  w*  =  9w/9x,  etc. 
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PRODUCTION  OF  PRECISION  FOAM  PLASTIC  PRODUCTS 
BY  PRESSURE  MOLDING  OF  FOAM  SLABS  IN  A  HOT  MOLD 

D.A.  Baker 

Flight  Research  Laboratory 
National  Aeronautical  Establishment 


SUMMARY 


This  report  descries  a  method  of  molding  model  Crash  Position  Indicators, 
which  will  bo  used  in  supersonic  wind  tunnel  tests,  from  blocks  of  styrofoam.  It  also 
briefly  describes  methods  that  were  tried  before  reaching  a  satisfactory  one,  and  lists 
a  few  arsas  where  this  type  of  manufacture  might  prove  valuable. 


INTRODUCTION 

In  testing  the  Crash  Position  Indicator  model  aerofoils  In  the  supersonic  wind 
tunnel  it  was  necessary  to  develop  a  method  of  manufacturing  a  lightweight  aerofoil  that 
could  be  made  quickly,  be  tough  enough  to  withstand  the  forces  from  supersonic  wind 
tunnel  speeds,  and  yet  be  light  enough  so  that  weight  could  be  added  to  bring  the  centre 
of  gravity  to  any  required  position  and  still  keep  the  density  of  the  aerofoil  within  the 
required  region. 


MOLDING  METHODS  AND  MATERIALS 

Five  different  molding  methods  and  materials  were  tried.  These  are  listed 

below, 

"Pour  In  Place"  Urethane  Foam 

This  material  formed  a  lightweight  model  but  the  edges  of  the  model  were 
too  delicate  and  friable.  It  was  difficult  to  remove  the  light  density  units  from  the  mold, 
and  proper  venting  of  the  mold  was  difficult. 

Syntactic  Foam 

This  material  is  made  from  microscopic  hollow  spheres  of  phenolic  resin 
with  a  polyester  resin  as  the  binding  material.  It  formed  a  good  strong  model  in  the 
deneity  range  of  15  to  20  pounds  per  cubic  foot. 

Polyethelene  Molding  Powders 

Models  were  formed  using  polyethelene  molding  powder  and  the  rotational 
molding  method.  This  provided  a  hollow  aerofoil  with  a  thin  tough  surface,  but  it  was 
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necessai’y  to  fill  the  cavity  with  foam  in  order  to  make  it  rigid  enough  for  supersonic 
speeds.  The  total  weight  was  prohibitive  and  the  method  dropped. 

Expandable  Polystyrene  Beads  (Dylite) 

This  material  will  expand  and  bond  together  when  heated.  It  provided  a 
lightweight  model  that  could  be  kept  well  within  the  weight  limits,  but  the  edges  of  the 
beacon  were  quite  weak. 

Polystyrene  Foam  (Styrofoam) 

Models  were  formed  by  pressure-molding  blocks  of  styrofoam  in  a  heated 
mold.  This  method  proved  to  be  a  successful  way  of  manufacturing  a  lightweight  aero¬ 
foil  with  tough  edges.  The  density  of  the  aerofoil  was  easily  controlled  by  controlling 
the  size  and  density  of  the  block  of  styrofoam  used  in  the  molding  process.  The  density 
of  the  edges  could  be  controlled  by  bevelling  the  edges  of  the  block  of  foam  (Fig.  1). 


METHOD  USED  TO  ’’HEAT  FORM"  POLYSTYRENE  FOAM 

The  recommended  method  to  "heat  form"  polystyrene  foam  is  to  soak  the 
foam  unit  at  a  temperature  slightly  below  the  heat  distortion  temperature  of  175°F, 
until  the  foam  block  is  a  uniform  temperature  throughout.  The  mold  is  heated  to  slight¬ 
ly  above  this  temperature,  to  190°  -  200°F.  The  heated  foam  block  is  then  placed  in 
the  mold  and  formed  to  shape  by  pressure. 

This  method  was  tried,  but  it  was  found  that  a  beacon  with  a  tougher  skin  and 
edge  was  obtained  if  the  foam  was  left  at  room  temperature  and  the  mold  temperature 
increased  to  220°  to  230°F  and  immediately  quenched  in  cold  water  when  the  foam  block 
was  molded  to  shape.  This  method  tended  to  compress  the  outer  surface  of  the  foam 
only,  and  produced  a  tough  skin  with  an  almost  solid  resin  edge. 


FOAM  BLOCKS 

The  foam  block,  for  a  square  beacon,  was  cut  as  shown  in  Figure  1.  The 
size  of  the  block  of  foam  was  made  larger  than  the  dimensions  of  the  mold.  Length  and 
breadth  were  increased  by  1/8  inch  in  order  to  take  care  of  misalignment  of  the  foam 
block  in  the  mold.  The  height  of  the  foam  block  varied  with  the  density  required  in  the 
finished  model.  It  was  found  that  the  minimum  height  of  foam  had  to  be  at  least  1/16  inch 
greater  than  the  depth  of  the  mold  to  insure  complete  molding  contact. 

Two  densities,  1.  8  and  4.3  pounds  per  cubic  foot,  were  tried  in  the  molding 
process.  The  height  of  the  block  of  foam  varied  from  1/2  to  1-1/4  inches.  The  densi¬ 
ties  of  the  aerofoil  ranged  from  3.86  to  21.5  pcf  (Fig.  2). 


MOLDS 

Two  types  of  molds  were  used  in  the  production  of  the  model  aerofoils.  The 
first  type,  Figure  3,  was  cast  from  Devcon  C,  which  is  a  mixture  of  80%  aluminum 
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powder  and  20f7c  epoxy  resin.  This  is  a  high  temperature  molding  compound  that  will 
withstand  5000  psi  at  400°  F,  with  a  shrinkage  on  curing  of  0.0001  inch  per  inch. 
(Manufacturer's  claim.) 

The  second  type  of  mold  was  carved  from  solid  aluminum.  It  had  a  better 
heat  transfer  than  the  plastic  aluminum  compound,  and  therefore  provided  a  quicker 
production  method.  This  type  of  mold  was  used  on  the  later  models.  The  model  pro¬ 
duced  from  the  above  mold  is  shown  in  Figure  4. 


CONCLUSION 

This  method  of  molding  provides  a  unique  way  of  controlling  the  density 
throughout  a  molded  product.  The  block  of  foam  that  is  to  be  pressure-molded  can  be 
roughly  shaped  to  produce  required  densities  in  different  areas  of  the  molded  unit.  The 
surface  density  can  be  controlled  by  the  temperature  of  the  mold.  A  thin  layer  of  solid 
resin  can  be  created  on  the  surface,  or  it  can  remain  a  high-density  foam  surface, 
depending  on  the  temperature  of  the  mold.  The  sharp  edges  of  the  molded  unit  can  be 
compressed  to  an  almost  solid  resin  or  high-density  foam  depending  on  the  thickness, 
at  the  edge,  of  the  foam  block  that  is  being  heat  and  pressure  formed. 

It  lends  itself  particularly  to  the  molding  of  any  section  that  requires  a  strong 
thin  edge,  such  as  the  trailing  edge  of  streamlined  struts,  wing  sections,  tail  sections, 
and  elevators  for  model  aircraft;  or  other  units  where  light  weight,  strength,  and 
toughness  are  requirements. 

This  method  can  also  be  used  to  mold  lightweight  filler  sections  for  fibreglas 
propeller  blades,  ailerons,  tail  sections,  etc. ,  for  full-scale  aircraft. 

It  provides  a  void-free  homogeneous  foam  unit  with  a  solid  resin  or  high- 
density  foam  skin  from  one  molding  operation. 


FOR  2"x  2  BEACON 
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FIG.  2*  MODEL  DENSITY  vs  SIZE  OF  FOAM  BLOCK  CHART 
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NAE  FLIGHT  IMPACT  SIMULA  TOR 


J.W,  Noonan  and  J.  B.  R.  Heath 
Structures  and  Materials  Laboratory 
National  Aeronautical  Establishment 


1.0  INTRODUCTION 

As  aircraft  become  larger,  fly  faster,  and  continue  to  extend  their  field  of 
operation,  so  too  must  test  and  development  facilities  match  this  extension  af  capability 
to  maintain  a  high  degree  of  safety  in  operation.  Of  all  the  hazards  to  flight  safety  the 
one  categorized  as  "bird  hazards  to  aircraft"  has  captured  the  imagination  of  almost 
everyone  involved  with  aircraft  and  is  currently  receiving  close  and  detailed  attention 
by  air  regulatory  bodies,  airline  operators,  and  military  air  commands.  The  progress 
reported  at  the  recent  World  Conference  on  Bird  Hazards  to  Aircraft  held  in  Kingston, 
Ontario,  in  September,  19G9,  gave  evidence  that  while  a  significant  decrease  in  bird 
incidents  involving  aircraft  in  the  vicinity  of  airports  had  been  achieved,  the  incidence 
of  en  route  encounters  remained  virtually  unaffected  and,  in  fact,  showed  signs  of  being 
on  the  increase.  It  would  seem  evident  that  the  only  measures  offering  any  real  assur¬ 
ance  of  reducing  the  en  route  hazard  arc  those  of  "bird-proofing"  of  the  aircraft. 

Bird-proofing  the  aircraft  consists  in  increasing  the  strength  and  the  energy- 
absorbing  capabilities  of  those  vulnerable  parts  of  the  aircraft  that  arc  likely  to  be  sub¬ 
jected  to  bird  impacts,  to  a  point  where  the  statistical  probability  of  a  catastrophic 
failure  owing  to  a  bird  strike  is  reduced  to  that  which  experience  has  dictated  as  an 
acceptable  level.  The  areas  considered  to  be  vulnerable  are  the  windshield^,  the  lead¬ 
ing  edges  of  the  empennage  surfaces,  and  openings  where  ingestion  into  the  engines  can 
occur;  wing  structures  are  generally  considered  to  have  sufficient  depth  to  sustain  a 
bird  strike  without  suffering  serious  damage.  A  great  deal  of  useful  design  information 
can  be  generated  by  simulating  bird  impacts  on  various  targets  over  an  appropriate 
range  of  velocities  and  environmental  conditions;  force-time  relations  can  be  derived 
for  a  range  of  bird  sizes;  failure  modes  can  be  identified  for  different  materials  and 
combinations  of  materials  and  the  ability  to  withstand  penetration  by  bird  sizes;  and 
velocities  spelled  out  in  airworthiness  regulations,  where  they  exist,  can  be  demon¬ 
strated. 


The  addition  to  the  facilities  of  NAE  of  the  Flight  Impact  Simulator,  popularly 
known  as  the  "bird  gun",  has  provided  a  means  for  simulating  bird  strikes  that  is  useful 
in  providing  information  applicable  to  designing  "bird-proof"  aircraft  and  for  determin¬ 
ing  the  penetration  velocity  for  critical  locations  on  existing  aircraft. 


11.0  THE  CH'N 

The  core  of  the  simulator  is  the  pneumatic  cannon  or  compressed  air  gun  and, 
as  it  has  been  described  in  the  NAE  Report  I.TR-ST  281  from  aspects  of  the  operating 
principles  and  detailed  design,  a  brief  description  here  should  be  sufficient.  The  gun, 
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a  photograph  of  which  is  shown  in  Figure  1  and  a  schematic  diagram  in  Figure  2,  is 
approximately  70  feet  long  and  consists  of  a  60-foot3  reservoir,  a  transition  section,  a 
pressure  step  chamber  that  also  functions  as  the  breech,  and  a  40-foot  barrel  10  inches 
in  diameter;  attached  to  the  muzzle  end  is  a  sabot  catcher.  The  supporting  frame  is 
pivoted  on  a  horizontal  axis  at  the  reservoir  end  to  permit  elevation  of  the  barrel  to 
facilitate  aiming  at  a  target  (Fig.  3). 


3.0  THE  SITE 

The  simulator  is  located  on  the  grounds  of  the  National  Aeronautical  Estab¬ 
lishment  at  Uplands  Airport  and,  as  can  be  seen  in  the  aerial  photograph  shown  in  Fig¬ 
ure  4,  it  is  housed  in  the  end  unit  of  a  row  of  Butler-type  buildings  that  lie  in  a  northerly- 
southerly  direction.  At  the  north  end  of  the  building  a  concrete  slab,  50  feet  x  75  feet, 
in  which  steel  tie-down  rails  are  embedded  forms  the  target  area;  a  removable  weather 
shelter  protects  the  central  part  of  the  area.  Approximately  150  feet  from  the  end  of 
the  building  a  semi-circular  earthen  butt  has  been  formed  to  contain  ftying  debris  result¬ 
ing  from  the  gun  operations.  As  can  be  seen,  the  target  area  is  readily  accessible  by 
service  roads  and  can  accommodate  complete  aircraft  when  the  weather  shelter  is 
removed. 


4.0  OPERATIONAL  DETAILS 

4 . 1  The  Sabot 

The  10-inch  gun  barrel  permits  a  range  of  bird  sizes  to  be  used  in  the  simu¬ 
lation  of  strikes.  International  airworthiness  codes  have  established  a  4-pound  bird  as 
a  standard  weight  for  impact  tests;  consequently  the  gun  is  extensively  used  for  this 
bird  weight  It  is  necessary  to  package  the  bird  in  a  "sabot”  that  fits  the  barrel  and  will 
contain  the  bird  without  breaking  up  until  it  leaves  the  muzzle  of  the  gun.  P’igurc  5  shows 
the  component  parts  of  the  sabot  for  a  4-pound  bird.  Because  of  the  relatively  large 
bore  of  10  inches,  the  present  sabot  is  only  usable  to  a  charge  pressure  of  50  psi;  this 
produces  a  velocity  of  420  knots  and  caters  for  almost  all  of  the  existing  requirements. 

4. 2  Bird  Types 

Both  synthetic  and  real  birds  are  used  in  the  simulator.  The  synthetic  bird 
is  made  from  a  stable  gel  formulation  consisting  of  cellulose  gum  (sodium  carboxy 
methylcellulose),  aluminum  acetate,  and  98 7i  water.  The  gel  is  cast  in  a  cylindrical 
form  about  5  inches  in  diameter  and  9  inches  long,  the  approximate  dimensions  of  the 
real  bird  package,  and  packaged  in  a  plastic  envelope  enclosed  in  a  cotton  bag.  Fig¬ 
ure  6  shows  a  real  and  a  synthetic  bird  ready  for  packaging;  no  attempt  has  been  made 
as  yet  to  simulate  any  part  of  the  bird's  skeleton  in  the  gel  bird. 

When  real  birds  are  used,  they  are  selected  from  domestic  fowl,  electrocuted, 
packaged  in  a  cotton  bag  and  stored  in  a  deep  freeze  unit  until  required  for  use.  The 
birds  are  withdrawal  24  hours  before  use  to  allow  them  to  reach  room  temperature. 

4.3  Diaphragms 

After  an  initial  assessment  of  scribed  2S-0  aluminum  diaphragms  and  limited 
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testing  of  0.0035-inch  aluminum  foil  in  single  and  multiple  sheets,  their  use  was  dis¬ 
continued  in  favour  of  the  polyester  film  known  as  mylar.  Diaphragms  cut  from  this 
material  proved  to  give  sharp,  repeatable  rupture  values  over  the  fairly  wide  range  of 
thicknesses  in  which  the  material  is  available.  Figure  7  gives  a  table  of  rupture  values 
for  the  10-meh  diaphragms,  and  Figure  8  illustrates  short-term  creep  characteristics 
of  the  sheet  material  at  room  temperature  for  the  10-inch  diameter  diaphragm;  it  can 
be  seen  that  the  material  can  be  held  under  pressure  for  a  considerable  time  without  the 
danger  of  an  unscheduled  firing  taking  place. 

4.4  Gun  Firing  Sequence 

After  inserting  the  bird  package  and  installing  the  diaphragms,  the  reservoir 
and  step  chamber  are  pressurized  in  the  ratio  of  2:1  to  the  required  charge  pressure, 
each  diaphragm  supporting  1/2  the  reservoir  pressure;  depressing  the  operating  button 
energizes  a  solenoid  valve,  allowing  the  pressure  to  bleed  off  from  the  step  chamber, 
and  the  diaphragms  are  ruptured  successively.  The  stored  energy  of  the  compressed 
air  in  the  reservoir  accelerates  the  bird  package  down  the  barrel,  the  sabot  is  arrested 
by  the  catcher  at  the  muzzle,  and  only  the  bird  travels  to  the  target. 


5.0  INSTRUMENTATION 
5. 1  Velocity  Measurement 

The  velocity  of  the  bird  at  impact  is  the  single  most  important  parameter  to 
be  measured  and  considerable  care  needs  to  be  exercised  in  setting  up  a  velocity  meas¬ 
urement  system.  Prior  to  selecting  a  particular  method,  a  number  of  different  systems 
were  tried.  They  were: 

(i)  breaking  of  stretched  wires  set  in  the  path  of  the  bird  to  interrupt 
electrical  circuits 

(ii)  breaking  of  tapes  set  in  the  path  of  the  bird  causing  light  beams  to 
impinge  on  photoelectric  units 

(iii)  interruption  of  a  laser  beam  focused  by  an  optical  system  on  photo¬ 
electric  units 

(iv)  an  optical  system  using  a  relatively  wide  and  narrow  light  beam  focused 
on  photodiode  units  in  conjunction  with  electronic  counter-timer  units. 

Methods  (i)  and  (ii)  suffered  from  the  same  disadvantages,  in  that  wires  or 
tapes  had  to  be  installed  before  every  firing.  The  wires  or  tapes  had  to  be  stretched  to 
consistent  values  to  achieve  producible  results  and  this  was  difficult  to  accomplish 
over  the  temperature  range  of  operations  from  -25°F  to  +100°F.  In  addition,  the  hold¬ 
ing  frames  had  to  be  securely  anchored  to  prevent  movement. 

Method  (iii),  in  which  a  laser  beam  was  focused  through  a  beam  splitter  and 
a  penta  prism  to  produce  two  parallel  beams  normal  to  the  bird  path  (both  beams  being 
focused  on  photoelectric  cells  to  initiate  and  terminate  a  count  of  the  time  interval  when 
interrupted  by  the  passage  of  the  bird)  was  undoubtedly  a  singularly  accurate  one,  but 
was  found  to  be  unworkable  from  a  practical  standpoint.  Problems  of  maintaining 
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alignment,  and  the  extreme  sensitivity  of  the  device  by  which  a  false  count  could  be  gen¬ 
erated  by  something  as  small  as  a  raindrop,  caused  (iii)  to  be  ruled  out. 

Method  (iv)  sounds  complicated  but  is  in  actual  fact  extremely  simple.  It 
consists  of  an  ordinary  automotive-type  sealed  boam  spotlight  with  a  pnrnbolic  reflector 
that  is  housed  in  a  box  with  the  front  face  masked  off  to  produce  a  wide  thin  beam,  pro¬ 
ject'  J  normally  through  the  bird's  path,  to  focus  on  u  2-inch  collector  lens  that,  in 
turn,  is  focused  on  a  photodiode  unit  incorporated  in  an  electronic  trigger  circuit.  'Pwo 
units  are  installed  approximately  7  feet  apart,  and  interruption  of  the  light  beams  by  the 
bird  generates  signals  to  operate  a  digitized  electronic  counter-timer  to  record  the  time 
elapsed  l>ctwocn  signals.  When  the  count  is  terminated,  the  digital  display  of  the  time 
interval  is  locked  in  and  ennnot  be  affected  by  further  events  until  the  counter  is  reset 
to  zer  o.  Two  completely  separate  systems  arc  used,  ono  in  the  horizontal  and  the  other 
in  the  vertical  plane.  This  redundancy  was  provided  in  case  of  a  system  malfunction. 
The  outstanding  feature  of  this  system  is  that  approximately  50%  of  the  light  reaching 
the  collector  lens  must  be  interrupted  before  triggering  a  count  signal,  thus  minimizing 
the  chance  of  a  false  count  by  small  stray  objects;  at  the  50%  obscuration  level  a  signal 
is  generated  by  Interrupting  the  beam  within  a  band  of  0,005  Inch  In  width.  Other 
advantages  are  that  nothing  needs  to  be  installed  prior  to  cnch  firing,  nnd  since  the  bird 
produces  no  reaction  with  the  light  beam,  the  frame  holding  the  system  has  only  to  be 
lightly  restrained,  In  15  months  of  operation  no  failure  of  the  system  has  occurred  and 
maintenance  has  been  confined  to  the  occasional  replace’  tent  of  one  of  the  spotlights, 
Figure  9  is  a  graphical  representation  of  a  unit  of  the  system  and  Figure  10  shows  the 
complete  setup, 

5.2  High  Speed  Photography 

A  16-mm  high  speed  camera  of  the  revolving  prism  type  is  part  of  the  normal 
instrumentation  of  the  simulator.  The  maximum  speed  of  the  camera  is  7000  frames 
per  second,  which,  when  projected  at  normal  film  speeds,  slows  down  the  action  ap¬ 
proximately  300  times.  In  normal  practice  the  bird  impnet  is  photographed  at  4000 
frames  per  second,  as  increasing  the  film  speed  beyond  this  value  requires  n  very  sub¬ 
stantial  increase  in  the  already  high  level  of  Illumination.  An  electronic  unit  is  used 
to  control  the  camera  sturt  and  its  speed  to  insure  the  capture  of  the  event  on  the  film. 
Initiating  a  typical  firing  sequence  by  depressing  the  operating  button  causes  (i)  the 
camera  to  start  about  0.5  second  before  the  bird  reaches  the  target,  and  (ii)  the  sole¬ 
noid  to  bleed  off  the  step  chamber  pressure  approximately  0,2  second  after  initiation. 
Since  in  this  particular  sequence  it  takes  0.2  second  from  the  time  the  dinphrngm  is 
ruptured  and  the  bird  impacts  the  target,  the  event  is  recorded  when  the  camera  has 
reached  its  predetermined  speed.  Study  of  the  filmed  events  arc  of  great  value  in  as¬ 
sessing  the  behaviour  of  the  target  materials  and  the  break-up  of  the  bird  at  impact. 
Examples  of  high  speed  film  coverage  of  (a)  a  4-pound  bird  travelling  at  317  knots  in 
impact  with  the  centre  pnncl  of  a  CF100  bullet-proof  windshield  arc  s^own  in  Figure  11, 
and  (b)  an  8 -pound  bird  travelling  at  315  knots  in  impact  with  the  leading  edge  of  a  DCS 
vertical  stabilizer,  in  Figure  12. 


G.O  CALIBRATION 

It  is  necessary  to  be  able  to  predict  the  velocity  of  the  bird  before  firing;  the 
basis  for  this  prediction  is  the  relation  established  between  the  reservoir  charge  pres¬ 
sure  and  the  velocity  achieved  for  a  given  bird  weight  us  measured  by  the  velocity 
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measurement  system.  A  series  of  test  firings  of  4-pound  birds,  up  to  a  velocity  of 
700  it/sec,  was  carried  out.  The  air  pressure  in  the  reservoir  was  measured  by  a 
mercury  manometer  to  within  0.05  in  Hg  (0.025  psi);  the  total  package  of  sabot  and 
bird,  weighing  approximately  7.0  pounds,  was  accelerated  the  length  of  the  barrel,  and 
after  the  sabot  was  arrested,  the  bird  continued  its  flight  through  the  measuring  system, 
the  centre  of  which  was  approximately  9  feet  from  the  muzzle  of  the  gun.  This  was  the 
point  at  which  the  average  velocity  between  the  7-foot  spacing  of  the  photodiode  units 
was  recorded  Synthetic  gel  birds  and  real  bird  carcasses  were  both  used  and  no  dif¬ 
ference  in  velocity  (that  could  be  attributed  to  the  difference  in  bird  types)  was  estab¬ 
lished.  Figure  13  gives  the  relation  determined  between  reservoir  pressure  and  bird 
velocity,  as  carried  out  under  the  conditions  noted  above,  up  to  the  limit  of  the  capa¬ 
bility  of  the  present  sabot. 


T.o  DISCI  SSION  OF  PERFORMANCE 

The  theoretical  basis  for  the  design  of  the  gun  was  taken  from  a  simple  per¬ 
formance  theory  (Ref.  1)  in  which  it  is  assumed  that  air  stored  in  a  reservoir  of  volume 
V,  at  pressure  P,  propels  a  bird  package  of  weight  W  along  the  barrel  of  length  L  and 
area  A  by  adiabatic  expansion.  From  the  equation  of  motion  of  the  bird  package  of 
weight  W  when  it  has  travelled  a  distance  s  along  the  barrel  and  the  pressure  has 
reached  a  value  P2  with  an  opposing  pressure  P0  on  the  other  side  of  the  bird,  the  veloc¬ 
ity  u  of  the  bird  is 
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where  y  ratio  of  specific  heats  for  air.  By  assuming  that  P0  (atmospheric  pressure) 
remains  constant,  and  neglecting  leakage  and  friction  losses  from  (1)  and  (2) 
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where  Ve  is  taken  as  the  velocity  after  travelling  distance  L,  i.e.  the  theoretical  muz¬ 
zle  velocity.  Integrating  gives 
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where 


(S) 


and 


reservoir  volume  ^  i 
gun  barrel  volume  AL 


Taking  V0  as  the  actual  velocity  attained,  then  V0  =  K2V'C  where  K2  is  a  correction 
factor  that  varies  with  velocity  and  the  particular  gun  design.  For  the  NAE  gun,  the 
values  of  Va  obtained  experimentally  allowed  K2to  be  determined,  assuming  negligible 
difference  between  muzzle  velocity  and  the  velocity  as  measured  9  feet  from  the  muz¬ 
zle.  Figure  14  shows  values  of  K2  plotted  against  V0  over  the  range  of  velocities  ai 
which  the  gun  has  been  operated;  it  is  evident  that  the  gun  becomes  less  efficient  at 
higher  velocities.;  Extrapolation  of  the  calibration  curve  shown  in  Figure  13  gives  on 
indication  that  the  maximum  velocity  achievable  at  a  charge  of  100  psig  would  be 
900  ft/sec  for  the  standard  4-pound  bird  as  fired  in  the  present  sabot.  The  provision 
made  for  evacuating  the  air  in  the  gun  barrel  should  make  it  possible  to  achieve  a  VQ  of 
1000  ft/sec  or  higher.  However,  the  present  sabot  is  limited  to  a  velocity  of  700  ft/sec 
above  which  it  is  not  completely  arrested  by  the  catcher  and  the  possibility  exists  that 
portions  of  the  sabot  could  reach  the  target.  New  sabot  designs  with  greater  energy¬ 
absorbing  characteristics  are  under  consideration,  as  are  other  means  of  dissipating 
the  energy  of  the  sabot  to  avoid  contacting  the  target. 


A  significant  amount  of  the  energy  of  the  air  stored  in  the  reservoir  under 
pressure  is  expended  in  accelerating  the  3-pound  sabot  in  order  to  achieve  the  required 
bird  velocity.  It  is  interesting  to  note,  however,  that  to  minimize  the  sabot  problem  by 
reducing  the  barrel  diameter  to  a  smaller  size  and  using  the  4-pound  bird  without  the 
sabot  results  in  a  lower  predicted  velocity  for  a  given  reservoir  pressure.  Consider¬ 
ing  the  NAE  gun  where 

W  (bird  +  sabot)  =  7.0  lb 
L  (barrel  length)  =  40  ft 

A  (barrel  area)  =  78.54  in2 

V,  (reservoir  volume)  =  60  ft3 
P,  (reservoir  pressure)  =  45  psia 
P0  (atmospheric  pressure)  =  15  psia 


from  (7) 


from  (6) 


2,75 


C,  =  0.804 


■w 
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and  from  ;'5j 


Vc 


783  ft/s<;c 


Reducing  the  barrel  diameter  to  6  inches  and  the  weight  of  the  bird  package  to  4  pounds, 
i.e, 

VV  *  4.0  ib 

A  =  28,27  in1 
V, 

K  -  — -  7.65 


and 


C,  -  0.921  I 

Vc  --  694  ft/sec  j 


At  45  psia  (30  psig)  the  actual  velocity  achieved  with  the  10 -inch  diameter  barrel 
and  vhe  7-pound  bird  package  war  560  ft/sec  and  gave  a  value  for  K2  of  0.715,  which 
represents  an  increase  over  values  for  6-  inch  diameter  barrel  configurations  reported 
for  ether  installations  (Ref.  1,  3>. 

The  operation  ox  the  simulator  in  conjunction  with  the  current  sabot  design 
has  proved  to  be  accurate  and  efficient  in  the  velocity  range  from  250  to  700  ft/sec,  the 
range  o.(  most  interest  tor  bird  impact  simulation,  and  serves  the  purpose  for  which  it 
was  designed.  Further  development  to  extend  the  velocity  to  the  1000-ft/sec  level  will 
be  undertaken  in  the  order  in  which  priorities  are  established. 


8,0  TYPICAL  IMPACT  SIMULATIONS 


Figure  15  shows  the  results  of  impacts  on  the  leading  edge  of  a  DC8  horizon¬ 
tal  stabilizer.  The  left-hand  impact  was  made  by  a  4-pound  gel  bird  travelling  at 
380  knots  and  the  right-hand  by  a  real  bird  carcass  travelling  at  357  knots.  Superficial¬ 
ly  the  two  impacts  look  somewhat  the  same  but,  v.'hereas  the  gc!  bird  merely  penetrated 
the  leading  edge  skin,  the  real  bird  penetrated  the  leading  edge  skin,  the  main  spar  web, 
and  some  of  the  denser  portions  of  the  bird  penetrated  the  rear  spar;  the  fact  that  the 
gei  bird  was  travelling  at  a  greater  velocity  than  the  real  bird  indicates  the  magnitude 
of  the  problem  of  producing  an  acceptable  synthetic  substitute  for  the  real  bird. 


Figure  16  shows  a  closeup  view  of  the  results  of  the  8-pound  bird  impact  with 
the  leading  edge  cf  a  DC8  vertical  stabilizer,  shown  at  high  speed  in  Figure  12.  The 
bird  was  travelling  at  315  knots  and  penetrated  the  triple  leading  edge  skins,  the  main 
spar  web,  several  canted  ribs,  and  substantial  portions  of  the  bird  contacted  the  rear 
spar,,  Ibis  particular  shot  was  made  in  response  to  the  FAA  invitation  to  comment  on 
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their  proposal  that  the  leading  edges  of  the  empennage  surfaces  of  transport  aircraft 
should  be  required  to  demonstrate  the  ability  to  sustain  an  a  -pound  bird  impact  at  cruise 
velocity  at  sea  level  (Ref..  2);  comments  have  been  forwarded  to  FAA.  Figure  17  shows 
the  results  of  the  4-pound  bird  in  impact  with  the  centre  panel  of  a  CF100  windshield, 
shown  at  high  speed  in  Figure  11.  The-  multi-layered  glass  and  plastic  p::nel  were  able 
to  withstand  penetration  by  the  bird,  but  the  east  magnesium  frame  was  broken  in  three 
places;  the  most  prominent  break  can  be  seen  at  the  top  left-1.,:  aJ  sub  of  the  frame. 

The  fifth  frame  from  the  top  in  the  high  speed  p.iotograuhs  shown  in  figure  11,  is  also 
worthy  of  note.  It  shows  that  the  liquid-like  flow  of  the  bird  on  contact  can  cause  high 
hydraulic  forces  to  develop,  and  in  this  case  bird  fluids  and  tissue  extruded  lx* tween  the 
top  of  the  centre  panel  and  the  metal  frame,  even  extruding  through  torn  rivet  holes  on 
the  inside,  an  indication  of  an  extremely  important  detail  to  be  taken  into  consideration 
at  the  design  stage. 


9.0  CONCLUSIONS 

The  XAE  Flight  Impact  Simulator  has  been  developed  into  a  reliable  impact 
loading  tool  operating  through  a  velocity  range  from  250  to  700  ft,;  sec,  using  4-  and 
8-pound  birds  as  projectiles.  Performance  over  range  has  exceeded  the-  designed 
value  for  K2 of  0.58,  i.e,  the  actual  velocity  \'0  -  0.58  Vc  predicted  originally  has 
been  exceeded  by  approximately  12'  a  value  of  0.  70  has  been  achieved  for  K2  and 
VQ  -  0.70  Vc  at  700  ft/sec.  The  \ elocny  can  be  predicted  to  within  -3fr  and  measured 
to  w>thin  =0.5fo.  Further  development  to  increase  the  velocity  to  1000  ft/sec  will  be 
undertaken  as  allocation  of  priorities  permits. 
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MYLAR  DIAPHRAGMS 


SHORT  TERM  CREEP  CHARACTERISTICS  OF  0  003-INCH  MYLAR  DIAPHRAGM 


(a)  TRIP  BAND  SHOWN  IN  VERTICAL  PLANE 


(b)  MINIMUM  OF  50%  LIGHT  OBSCURED  TO  TRIGGER  SIGNAL 


(c)  LESS  THAN  50%  DOES  NOT  TRIGGER  SIGNAL 


FIG. 9:  VELOCITY  MEASURING  DEVICE 
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CURRENT  PROJECTS 


Much  of  the  work  in  progress  in  the  labora¬ 
tories  of  the  National  Aeronautical  Establishment  and 
the  Division  of  Mechanical  Engineering  includes  cali¬ 
brations,  routine  analyses  and  the  testing  of  proprie¬ 
tary  products;  in  addition,  a  substantial  volume  of 
the  work  is  devoted  to  applied  research  or  investi¬ 
gations  carried  out  under  contract  and  on  behalf  of 
private  industrial  companies . 

None  of  this  work  is  reported  in  the  follow¬ 


ing  pages . 
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ANALYSIS  LABORATORY 


AVAII.AIII.K  KACIUTIE8 


HiIh  laboratory  Hum  analysis  nnd  simulation  facilities  avnllnblo  on  nn  open-shop  basis,  Enquiries  are  especially  cncour- 
oHnrUv  Pr?  t'''lVnnt  m“y  ut,m.“  tl!t'  ,"cllltlc''  ln  a  novel  »nd/or  puitleularly  effective  manner,  Such  projects  are  given 
.MM  I  .  i.r'j'.  n  f  y  ?Up,J0rM,<?  W!th  BM,",t,ncc  frorn  InlHiratory  poraonnot.  The  facilities  arc  especially  suited  to  sys¬ 
tem  di sign  studies  and  scientific  data  processing.  Information  Is  available  upon  request,  1  y 


EQUIPMENT 

1,  An  Electronic  Associates  090  IIYItltll)  COMPUTE)!  consisting  of  the  following! 

(a)  EAl  040  digital  computer 

-  I0K  memory 

-  eard  reuder 

-  high  speed  printer 

-  disc 

-  digital  plotter 

(b)  EAI  OHO  analogue  enmpulnr 

-120  amplifiers 

-  non-linear  elements 

-  x-y  pen  recorders 

-  0-chunnel  strip  chart  recorder 

-  large  screen  oscilloscope 

(c)  EAl  093  Interface 

-  18  dlgltal-to-annlogue  converters 

-  32  analogue-to-dlgltal  converters 

-  Interrupts,  sense  lines,  control  linen 

2.  A  Precision  Instruments  Model  400  KM  tape  recorder.  HUG  atnndnrd  (10  Kllr.)  ns  well  ns  extended  bandwidth  IRIG 
ala  in  la  nl  (Dll  Kiln)  record  mill  iiliiylnirh  nlimii’iinlcsi  M  (lulu  ( ruphw  (1  Inch  tape)  plus  2  vnlro  ohnimols  or  7  rlnln  tracks 
ll/H  ll  llll  l«|!P);  Ift/ll)  |j||i|«i«  |MT  snioinl  III  I  I'll  Iflolinfi  P<  I  hl'IWl  l/||)fi  fl|)PPll: 


(IMNEIIA  I,  HTUDIEH 

Hindy  of  eomimllng  melhtxls  In  opllmal  control. 

Ilevelupteenl  of  slallsllenl  analysis  loelinli|ties  for  (lie  analysis  of  analogue  recorded  data,  Including  digital  and  hybrid 
Implementations, 

study  of  hybrid  computing  techniques  for  parameter  optimisation  In  dynnmle  systems, 

Study  of  the  problem  of  decoupling  multivariable  syslem*  by  slab’  vnrlnble  feedback. 


A PP  liICATIONH  STUDIES 

In  collaboration  with  the  Plight  Itesenrch  laboratory  of  NAE,  a  study  Is  being  made  of  magnetic  fields  generated  by 
surface  waves  and  Interval  waves  In  oceans. 

In  collaboration  with  the  Flight  Itesenrch  bnbornlory  of  NAE,  ntnllntlrnl  analysis  techniques  are  being  developed  for 
the  analysis  of  clear  air  turladence  records, 

In  collaboration  wllh  Graham  E,  Crate  I.td,  anti  Trans  Canada  Pipe  bines,  a  hybrid  computer  simulation  of  a  natural 
gas  pumping  station  la  Itelng  developed,  Numerous  control  schemes  for  operating  the  eompressors  In  parallel  are 
being  Investigated, 

Development  of  digital  programs  for  Ihe  analysis  of  data  obtained  try  NIICs  Hlosclenees  InlHiratory  In  Ihe  study  of  the 
effects  id  low  lenqiernture  on  rats, 

ln  collaboration  wllh  Ihe  Mechanical  Engineering  Dept,  of  Queens  University,  a  hybrid  technique  Is  being  developed  for 
the  solution  of  the  partial  differential  equallon  describing  (lie  gas  pressures  in  the  du.’ts  of  a  Iwo-stroke  engine. 

In  rollalsirnllon  with  Ihe  Engine  laboratory  of  IIME,  a  hybrid  computer  simulation  of  a  free-plston-acllvnled  hydraulic 
pump  Is  being  studied.  The  simulation  Is  isdng  used  lo  establish  ihe  engine  geometry  for  efficient  operation  over  wide 
hydraulic  power  levels. 

Development  of  programs  louse  a  magnetic  disk  ns  virtual  memory  for  the  EAl  M0  computer 

In  cullnliorntltm  wllh  Canadian  Aviation  Eleelronlea  In  Montreal,  dlgltlrlng  of  data  for  a  fully  automatic  compensation 
syslem  fo'  magnetic  anomaly  detection  Is  being  performed. 

Development  of  programs  lo  produce  a  parallel  perspective  drawing  of  a  function  of  two  variables,  using  the  digital 
plotter.  Also  some  study  on  the  elimination  of  hidden  lines  In  persperllve  drawings  and  on  contour  mapping  of  functions 
o(  two  variables, 


PRtCtlllhb  rmu.  vuUK 
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nwiT.ulT1  °U'  L'urv.°-flulnK  computer  program  to  allow  fitting  of  nnnlytle  functions  of  many  types  to  export- 

mental  data,  .he  package  provide*  on-line  oiierator/eomputor  dialogue  as  well  as  on-llno  plotting  of  results, 

rmB  f°r  P?rfor,mlnK  P,owor  apcct™'  density,  erosB-spoclral  density,  auto-  and  cross-correlation  analyses, 
Engineering  '  °Ur  L‘r  lransform  algorithm,  arc  being  developed  for  the  Hydraulic  .Section  of  the  Division  of  Mechanical 


in  collaboration  with  the  Canadian  Pacific  Hallway,  an  uppralsat  study  of  the  riding  properties  of  the  CPU  80-foot  con¬ 
tainer  cars  Is  being  conducted  on  the  EAI  690  hybrid  computer  system. 


APPLICATIONS  STUDIES  HY  OTHERS 

The  Radio  and  Electrical  Engineering  Division  Is  using  the  EAI  040  digital  computer,  dlgllnl-to-analogue  and  analogue- 
to-dlgttnl  converters  In  processing  data  for  the  study  of  polarization  of  solar  emission. 

The  Plight  Resenrch  group  of  the  NAE  Is  using  the  EAI  040  digital  computers,  dlgllal-lo-analogue,  annlogue-lo-dlgltul 
convertors  and  Cnlcomp  plotter  In  ohtninlng  power  spectra  of  MAD  (Magnetic  Anomaly  Detection)  signals. 


CONTROL  SYSTEMS  LABORATORY 


INDUSTRIAL  CONTROL  1’HOHLRMS 

Investigation  of  industrial  systems  applications  of  fluidic  circuits. 

In  collaboration  with  the  Department  of  Energy,  Mines  anti  Resources,  an  Investigation  of  the  process  dynamics  and 
control  characteristics  of  an  electric  arc  furnace  for  processing  Iron  ore. 

Dynamic  modelling  of  electric  arc  and  oxygen  steelmnklng  processes. 

Investigation  of  the  process  dynamics  and  control  characteristics  of  u  copper  converter, 


I A  HCiE  SYSTEMS  STUDIES 

Investigation  of  the  possible  Influence  of  fresh  water  outflow  on  climate. 

HUMAN  FACTORS  ENGINEERING 

A  general  program  of  research  arid  development  In  the  Human  Factors  Engineering  field  that  Includes  the  following! 

Investigation  of  the  control  characteristic",  of  the  human  operator  ami  the  basic  phenomena  underlying  tracking  perform¬ 
ance. 

Investigation  of  the  nature  of  sensory  Interaction  In  human  perceptual-motor  performance. 

Investigation  of  the  factors  Involved  In  the  presentation  and  processing  of  Information,  particularly  In  relation  to  simu¬ 
lator  design. 

IIIO  LOCI  1C  A  L  KNOIN  E  Elt  INC. 

A  general  program  of  research  and  development  In  the  biological  engineering  field  that  Includes  the  followings 
Investigation  of  the  Implementation  of  feedback  control  In  living  organisms. 

Investigation  of  data  transmission  processes,  with  particular  reference  to  nerve  conduction  characteristics. 

Investigation  of  auditory  methods  of  monitoring  electrnphyslologleal  signals  In  general  and  the  electrooncephnlognnh  In 
particular. 

Development  of  depth  probes  for  the  study  of  electrical  activity  In  the  deep  structures  of  the  human  brain. 

Development  of  sleroo-taxle  ami  other  a|iparatus  for  neurosurgical  procedures. 

Development  of  a  phase  memory  filter  for  electroencephalograph  studies. 

Investigation  of  musclo  control  by  ordered  electro-stimulation. 


PATTERN  RECOGNITION 

Investigation  of  the  fundamentals  of  pattern  recognition. 

Development  of  techniques  for  the  Identification  of  biological  cell  populations,  fingerprints  etc. 
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MUD  DISPERSAL  I1Y  MICROWAVE  RADIATION 


— 


ENGINE  LABORATORY 


KHKE  PISTON  ENGINES 

Study  of  tho  free  piston  engine  and  Its  applications.  Mathematical  study  of  the  dynamic  and  thermodynamic  processes  In 
the  engine,  using  computer  simulation  techniques,  with  a  view  to  elucidating  starting  and  control  phenomena.  Experi¬ 
mental  assessment  of  the  operation  and  performance  of  the  free  piston  engine  In  tho  laboratory.  Study  of  tho  Injection 
system  most  suitable  for  the  peculiar  requirements  of  the  free  piston  ongino,  with  possible  application  to  nil  direct  In¬ 
jection  engines,  The  computer  simulation  technique  has  recently  been  extomlod  to  tho  study  of  hydraulic  power  take-off 
from  n  free  plalon  engine  through  a  pneumatic  coupling. 


DUCTED  KAN  AERODYNAMICS 


Aertxlynnmle  performance  study  of  highly  loaded  ducted  fans,  with  particular  reference  to  Inlet  and  outlet  distortion 
phenomena  ns  encountered  typically  by  V'TOI,  aircraft.  Analytic  study  of  the  flow  through  a  short  compressor  that 
operates  with  non-uniform  Inflow  and  outflow  conditions.  Experimentally  assessed  performance  of  n  fan-in-wing  model 
In  a  closed  test  section  wind  tunnel  over  u  range  of  Inflow  ratios,  angles  of  attack,  and  with  a  variety  of  Intel  shapes. 


V/STOI.  NOISE  STUDIES 

Study  of  the  mechanism  of  the  generation  and  suppression  of  noise  produced  by  ducted  fnns  for  VTOI,  aircraft.  Identi¬ 
fication  of  the  noise  sources  and  relating  the  strength  of  the  sources  to  the  physical  parameters  of  tho  system. 

iN- I'T.ltlllT  THRUST  METER  KOI,  ,IKT  ENGINES 


An  Investigation  concerned  with  the  development  of  a  meter  capable  of  indicating  continuously  the  gross  thrust  of 
a  Jet  engine,  being  carried  out  In  conjunction  with  a  Canadian  electronics  firm. 


CENTRIFUGAL  COM  I’ll  ESI’, OHS 


Design  anil  performance  Investigation  of  centrifugal  compressors.  Including  study  of  flow  phenomena  In  overnlr.c  model 
Impellers.  Detailed  study  of  devices  for  stabilising  the  flow  In  the  channels  of  a  complete  but  simplified  Impeller,  Test 
facilities  have  been  expanded  to  permit  Improved  correlation  between  flow  patterns  and  performance 


AXIAL  COMPRESSORS 

Preliminary  analytical  and  experimental  studies  of  small  axial  compressors. 


IOCOMOTIVE  DIESEI.  ENGINE  PHOtU.EMS 


In  co-oporntlon  with  the  Canadian  National  Railway  and  the  Canadian  Pacific  Railway,  an  Investigation  of  locomotive 
•Reset  engine  problents,  including  those  arising  from  the  use  of  Athabasca  Tar  Sands  crude  oil  ns  fuel.  Studies  of  wear 
processes  In  the  engine  employing  e  specially  developed,  new  method  of  spectrogrnphle  sampling  of  cylinder  oil.  Inves¬ 
tigation  of  several  new  types  of  lubri-atlng  oil,  oil  filters,  different  kinds  of  cylinder  liners  etc.  with  respect  to  engine 
wear.  Exploratory  Invesllgallon  and  lesla  for  planning  a  maintenance  monitoring  system  on  board  diesel  locomotives  In 
service 


EOAMED-CIAY  MA  TERIA  US 


Investigation  <>f  novel  light-weight  foamed-rlny  building  materials,  with  respect  toehemleal  eomposltlcmnnd  physical  prop¬ 
erties,  Application  Is  being  made  for  patents  on  the  process,  nfler  which  It  will  lx*  offered  to  Canadian  Manufacturers 


ENGINE  LABORATORY 
ISiON  OF  MECHANICAL  ENGINEERING 


-  75  - 


FLIGHT  RESEARCH  LABORATORY 


DESIGN  AND  DEVELOPMENT  TESTING  OE  A  CHASII  POSITION  INDICATOR  FOR  HELICOPTERS 

Experiments  on  models,  supplemented  by  theoretical  analyses,  are  being  conducted  for  the  purpose  of  evolving  an  Im¬ 
proved  crash  position  Indicator  for  helicopters. 


AIR  HORNE  REMOTE  SENSING  OE  MAGNETIC  PHENOMENA 


Experimental  and  theoretical  studies  relating  to  the  further  development  and  use  of  magnetic  airborne  detection  equip¬ 
ment.  Equipment  under  development  is  installed  on  a  North  Star  aircraft,  which  is  used  as  a  flying  laboratory  and  for 
preliminary  surveys  requested  by  various  Government  Departments  and  other  agencies. 

A  Iteech  Queenair  aircraft  is  also  being  equipped  and  operated  for  the  Department  of  Energy,  Mines  and  Resources  as  a 
prototype  for  future  general  aeromagnctlc  surveys. 


AIRBORNE  REMOTE  SENSING  USING  INFRA-RED  TECHNIQUES 

The  same  North  Star  aircraft  is  being  used  to  investigate  the  potentialities  of  airborne  Infra-red  techniques  In  hydro- 
logical  surveying,  soil,  permafrost  studies,  and  other  possible  applications, 


INVESTIGATION  OE  E LYING  QUALITIES  AND  CONTROL  SYSTEM  REQUIREMENTS  APPLICABLE  TO  V/STOL  AIRCRAFT 

Airborne  simulation  techniques,  using  helicopters  equipped  to  provide  variable  stability  and  control  properties,  are 
being  employed  to  explore  the  effects  of  the  numerous  parameters  involved  and  to  produce  data  that  are  directly  appli¬ 
cable  for  design  purposes.  Specific  investigations  are  also  being  conducted  for  aircraft  manufacturing  firms  and  other 
agencies  Whenever  possible  direct  comparisons  are  made  between  results  obtained  using  the  helicopter  equipped  as  an 
airborne  simulator  ami  those  obLiined  on  actual  VTOL  and  STOL  aircraft.  A  Bell  47G3B1  helicopter  with  vnriable 
eharaeteristlcs  in  four  degrees  of  freedom  Is  used  for  current  research,  A  Bell  205A1  helicopter  Is  being  developed  as 
an  airborne  simulator  with  variable  characteristics  in  six  degrees  of  freedom. 


INVESTIGATIONS  REIATING  TO  FOREST  EIRE  CONTROL  BY  AERIAL  METHODS 


Studies  of  various  factors  determining  the  effectiveness  of  aerial  fire  suppression  methods,  Including  theoretical  and 
experimental  work  on  the  behaviour  of  liquids  releas.  J  into  an  airstream,  operational  analyses,  and  the  Investigation  of 
aircraft  design  requirements.  Data  on  aircraft  behaviour  and  on  liquid  drop  patterns  are  being  obtained  from  flight 
experiments  with  a  number  of  typical  fire-bombing  aircraft. 


INVESTIGATION  OF  ATMOSPHERIC  TURBULENCE 


A  T-33  aircraft,  equipped  to  measure  wind  gust  velocities,  air  temperature,  wind  speed,  and  other  parameters  of 
interest  in  turbulence  research,  is  being  used  for  several  investigations.  These  include  measurements  at  very  low 
altitude,  in  clear  air  above  the  tropopruse,  in  the  neighbourhood  of  mountain  wave  activity,  and  near  storms.  Records 
are  obtained  on  magnetic  tape  to  facilitate  data  analysis.  Clear  air  turbulence  detection  methods  are  also  being  Investi¬ 
gated  The  aircraft  also  participates  in  ro-oporatlve  experiments  with  other  Canadian  and  foreign  research  agencies. 
A  second  T-33  aircraft  is  used  in  a  supporting  role  for  these  and  other  projects. 

A  USA F  RB-57F  aircraft  has  been  instrumented  ami  is  gathering  turbulence  data  on  routine  flights  above  50,000  feet. 


AIRCRAFT  OPERATIONS 


Various  studies  relating  to  aircraft  operations  are  made  from  time  io  time.  These  may  involve  such  matters  as  the 
provision  of  technical  assistance  during  accident  investigations,  the  analysis  of  particular  aspects  of  aircraft  behaviour 
in  operations,  or  the  preparation  of  recommendations  on  flight  recorder  requirements  and  specifications. 


INDUSTRIAL  ASSISTANCE 


Assistance  is  given  to  aircraft  manufacturers  and  other  companies  requiring  the  use  of  specialized  flight  test  equipment 
or  techniques. 


FLIGHT  INVESTIGATION  OF  SPRAY  DROPLET  RELEASE  FROM  AIRCRAFT 


A  Harvard  aircraft  Is  being  equipped  for  spray  nozzle  flight  experiments. 


FLIGHT  RESEARCH  LABORATORY 

ional  aeronautical  ESTABLISHMEN 


I 
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FUELS  AND  LUBRICANTS  LABORATORY 


iiydkcx;en-oxygen  engines 

Theoretical  studios  aimed  at  producing  the  host  design  features  for  hydrogen-oxygon  engines  in  multi-stage  vehicles. 

Investigation  of  heat  transfer  in  a  500-11)  thrust  rocket  combustion  chamber  using  a  water-cooled  chamber  burning 
hydrogen  and  oxygen. 

Testing  of  different  types  of  rocket  oropellant  injection  plates. 

Experiments  on  cryogenic  tankage. 


M  NDAMENTAt.  STl'DIES  OK  KKICTION,  I.l'BHICATION  AND  WEAK  PROCESSES 

Investigations  of  friction  and  wear  processes  Including  the  mechanism  of  adhesion  lx1  tween  non-conforming  metal  surfaces 
and  the  processes  Involved  with  transferred  flints  of  solid  lubricants. 

Theoretical  an* I  experimental  studies  of  a  modified  type  of  hydrodynamic  thrust  bearing  and  a  theoretical  study  of  a  novel 
form  of  gas  hearing. 


PRACTICAL  STL'DIES  ON  I.l'BRICATION,  FRICTION  AND  WEAR 

Assessment  of  wear  in  shotgun  barrels  with  shot  manufactured  from  different  materials. 

A  co-operative  program  for  the  assessment  of  Instrument  oils  and  lubricant  surface  coatings  in  the  bearings  of  miniature 
rotating  electrical  components. 


COM Bl'STION  RESEARCH 

Experiments  on  fuel  spray  evaporation 

EXTENSION  AND  DEVELOPMENT  OF  IAHORA TORY  EVALUATION 

Investigation  of  laboratory  engine  test  procedures  for  evaluation  of  oxidation,  dlspersancy,  and  thermal  stability 
characteristics  of  engine  oils. 

Development  of  laboratory  full-scale  axle  procedures  for  the  determination  of  the  anti-score  performance  of  hypoid 
gear  oils. 

PERFORMANCE  ASPECTS  OF  FEELS,  OII.S,  GREASES,  AND  BRAKE  FLL'ID 

Co-operative  investigation  covering  used  oil  analysis  anil  inspection  of  engines  from  Ottawa  Transportation  Commission 
buses  to  establish  realistic  oil  and  filter  change  periods. 

Evaluation  of  used  oils  from  railroad  diesels  to  establish  suitable  test  methods  and  condemning  limits. 

Development  of  a  specification  for  high  viscosity  index  hydraulic  oils  for  marine  use. 

Examination  and  evaluation  of  some  re-reflned  oils. 

Investigation  of  lalvoratorv  methods  for  predicting  flow  properties  of  engine  and  gear  oils  under  low  temperature  oper¬ 
ating  conditions. 

Investigation  of  laboratory  methods  for  predicting  low  temperature  flow  properties  of  diesel  and  heating  fuels  and 
assessment  of  their  suitability. 

Evaluation  of  methods  for  determining  undissolved  water  content  of  aviation  turbine  fuels, 
low  temperature  performance  of  hydraulic  oils  In  pump  systems. 

Investigation  of  lubricating  oil  performance  In  water-cooled  2-stroke  engines. 

Development  of  a  laboratory  method  for  evaluating  the  shear  stability  of  multigrade  motor  oils. 

Evaluation  of  additives  for  preventing  carburetor  Icing  In  light  aircraft  reciprocating  engines. 

Determination  of  fuel  system  Icing  Inhibitor  In  aviation  turbine  fuels, 

Performance  evaluations  on  domestic  furnace  fuel  oil  filters. 

Diesel  fuel  system  icing  studies. 

Investigation  of  the  Performance  of  synthetic  hypoid  gear  lubricants. 

M  ISC  E  LI  AN  MM'S  STl'DIES 

The  preparation  anil  cataloguing  of  Infra-red  spectra  of  compounds  related  to  fuels,  lubricants,  and  associated  products. 
The  application  of  Atomic  Absorption  spectroscopy  to  the  determination  of  metals  In  petroleum  products. 

Investigation  of  the  stability  of  highly  compressed  fuel  gases. 


500- lb  rocket  thrust  chamber  firing 


FUELS  AND  LUBRICANTS  LABORATORY 
DIVISION  OF  MECHANICAL  ENGINEERING 
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GAS  DYNAMICS  LABORATORY 


V/STOL  PROPULSION  SYSTEMS 

A  general  study  of  V/STOL  pro|Sjlslon  system  methods  with  particular  reference  to  requirements  of  economy  and 
safety. 

Investigation  of  a  vectored  lift/thrust  engine  arrangement  Involving  a  shrouded  fan  driven  by  a  partial  admission  turbine. 

Experimental  investigation  of  a  pod-mounted  VTOI.  fan  for  studies  Including  the  effects  of  flow  of  distortion  In  cross- 
flow  and  shroud  thrust  effects. 

Examination  of  V/STOL  propulsion  system  Intake  arrangements  Involving  lifting  surface  circulation  control. 
INDUSTRIAL  (IAS  TUH1HNES 

Investigation  of  gas  turbine  components,  particularly  for  Industrial  applications. 


INTERNAL  AERODYNAMICS  OK  DUCTS 

An  experimental  study  of  the  Internal  aerodynamics  of  ducts,  liends,  otxl  diffusers  with  particular  reference  to  the 
effect  of  entry  flow  distortion  In  geometries  Involving  changes  of  cross-sectional  area,  shape,  and  axial  direction. 


HEAT  TRANSFER  STUDIES 

The  study  of  heat  transfer  within  a  vertical  sealed  tube  (thermo-syphon)  in  which  working  fluid  is  boiled  In  the  lower 
section  ami  condensed  In  the  upper  section.  Several  practical  applications  for  this  device  are  being  studied. 


HYDROSTA  TIC  CAS  HEAR  INC. S 

Studies  of  hydrostatic  gas  bearings  to  develop  reliable  methods  of  predicting  hearing  performance  for  a  range  of  condi¬ 
tions  ami  configurations,  ami  to  evolve  suitable  teehnbiues  for  the  satisfactory  application  of  this  type  of  bcnrlng  In 
situations  where  the  special  properties  of  gas  liearlngs  recommend  their  use. 


ARC  PRODUCED  PIASMA  STUDIES 

A  general  Investigation  of  the  properties  and  flow  behaviour  of  thermally  Ionized  gases  produced  bv  arc  heating  on  a 
continuous  basis,  together  with  the  development  of  suitable  diagnostic  techniques  for  the  study  of  such  high  temperature 
gases. 

SHOCK  PRODUCED  PLASMA  STUDIES 

A  general  theoretical  and  experimental  investigation  of  the  production  of  high  temperature  plasma  by  means  of  shock 
waves  generated  by  electromagnetic  or  gasdynamic  means,  anil  the  development  of  diagnostic  techniques  suitable  for  n 
variety  of  shock  geometries  ami  the  study  of  physical  properties  of  such  plasmas. 

A  theoretical  and  experimental  study  of  strong  converging  ami  diverging  shock  waves,  produced  by  chemical  (explosive) 
ami  other  means,  ami  the  development  of  experimental  means  to  study  the  resulting  transient  plasma. 


HIGH  PRESSURE  LIQUID  JETS 

Water  jc'x  ranging  in  size  from  0.001"  to  0.015",  generated  by  pressures  in  the  range  of  10,000  to  100,000  psl,  arc 
capable  o.'  cutting  materials  such  as  paper,  cloth,  plastics,  wood,  masonry,  and  even  some  metals. 

Laboratory  work  Is  directed  towards  commercial  exploitation  by  various  industries  (c.g.  paper,  lumber,  leather,  gar¬ 
ments,  plastics,  etc.)  of  this  phenomenon,  and  to  the  detailed  study  of  the  phenomenon  itself  in  order  to  Improve  the 
efficiency  of  the  process  and  possibly  to  broaden  the  field  of  application. 

INDUSTRIAL  PROCESS  INSTRUMENTATION  AND  APPARATUS 

There  is  an  appreciable  effort,  on  a  continuing  basis,  directed  towards  industrial  assistance.  This  work  Is  of  an  ex¬ 
tremely  varied  nature  and,  in  general,  requires  the  special  facilities  and  capabilities  available  in  the  laboratory.  The 
following  examples  of  work  during  the  past  year  arc  typical:  altitude  performance  of  small  gas  turbine  engine;  hydraulic 
pump  tests;  cent-ifugal  compressor  tests  for  pipeline  service;  design  and  test  of  heaters  for  laying  continuous  rails 
on  Canadian  railways;  assistance  in  alr-bcarlng  design  for  Industrial  applications;  instrumentation  for  copper  and 
Iron  smelter  applications;  design  and  experimental  verification  of  pressure  loss  data  for  large  air  ducts  In  smelterst 
industrial  and  aero  turbine  development. 

Current  co-operative  projects  with  manufacturers  Include: 

(a)  Experimental  development  of  gas  turbine  combustion  chambers 

(b)  Theoretical  analysis,  design,  and  development  of  a  gas  bearing  spindle  for  the  textile  Industry 

(c)  Experimental  assessment  and  development  of  cushion  and  propulsion  schemes  for  track  air  cushion  vehicles. 
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H/GH  SPEED  AE RODYNAM/CS  LABORATORY 


CONTROL  OF  A  TL  RBI  LIST  BOl  SDAHY  LAVLK  IN  A  THRU  -DIMf  NsjDNAI.  SJhHK  WAVE  Blit  SDAHY  LAYER 
INTERACTION 

The  j-tn  •  5- in  bloinkmu  «i»i  tunnel  Is  t>«:ng  used  to  utni^tr  the  ibretwilir*  nsiueal  interaction  between  a  glancing, 
oblique  shook  live  and  a  turbulent  buundarv-Ltw-r  flow  along  <  tiat  wait  to  a  seeoud  phase,  the  boundary- layer  (law  in 
the  threv-di’nrosioaal  interaction  region  wilt  t*  rr-tncfeunl  tv  ungintiai  air  blowing. 


HYPERSONIC  INTERNAL  i'W  STID1ES 

The  fundamental  aspects  at  tcu-nul  flows  are  being  investigated,  to  the  theoretical  phase  of  this  prefect  the  InvtsUd 
flow  fields  of  intakes  with  arbitrary  lip  prattle  at  zero  incidence  arc  computed,  using  an  imploding  shock  analogy  and  a 
itwihi-sirip  shock  layer  tivorv .  Viscous  interaction  Ls  also  inci*nhd.  A  mixed  method  using  finite  differences  and  the 
.netbed  of  characteristics  wiil  be  developed  for  the  internal  flow  at  angle  of  incidence.  laminar  three-dimensional 
boundary  -layer  growth  will  also  tv  investigated.. 


SCHUERES  SYSTEM 

All  four  transonic  section  windows  have  been  installed  and,  using  a  matching  source  plate,  made  tc  correspond  to  the 
new  optical  path  length.  Good  photographs  have  been  takrn  through  the  transoaie  section  under  normal  testing  conditions, 
using  the  partially-completed  controller. 

Assembly  of  the  reground  supersonic  windows  is  proceeding,  a  film  identification  numbering  system  is  being  added,  and 
vari--_s  aspects  of  the  controller  are  being  completed  for  a  planned  final  Installation  to  February., 


NUMERICAL  SOLUTIONS  OF  LAMINAR  AND  Tt'RBl  LEST  BOUNDARY  I  .AYERS 

The  boundary- toyer  equations  for  a  two-dimensional  in<  omprcssible  flow  have  been  solved  by  a  finite  difference  method. 
For  the  turbulent  boundary  layer,  the  mixing  length  hypothesis  is  used  to  relate  the  turbulent  shear  stress  to  the  mean 
flow.  The  method  is  being  extended  to  include  the  solutions  of  boundary- layer  flows  to  a  compressible  fluid. 


NUMERICAL  SOLlTiONS  OF  THE  FLOW  FIELD  FOP  CONICAL  BODIES  D>  A  St  PERSONIC  STREAM 

A  numerical  procedure  for  solving  the  prob  cm  of  steady  supersonic  inviscid  flow  around  smooth  conical  bodies  h_j  been 
developed.  The  procedure  solves  the  elliptic  partial  differential  equation  that  define  the  conical  flow  between  the  body 
and  the  shock.  Results  have  been  obtained  for  circular  cones  up  to  relative  incidences  of  about  1.4,  including  some 
cases  for  incidences  beyond  a  critical  value  i  which  the  entropy  singularity  is  off  the  body  surface.  Also,  results  have 
been  obtained  for  cooical  bodies  with  non-circular  cross  sections  such  as  elliptic  copes.  Results  for  circular  cones  are 
presented  in  Agzrdograph  No.  137,  Sente  mb'."-  1969, 


TWO-DIMENSIONAL  TRANSONIC  TEST  FACILITY 

This  15-in  x  60-in  facility,  briefly  described  on  a  "would  be"  basis  in  previous  issues  of  this  bulletin,  was  installed 
and  operated  for  the  first  time  in  October.  An  ambitious  calibration  program  was  originally  planned,  but,  for  various 
reasons,  this  ha-  to  be  reduced  to  a  bare  minimum.  However,  enough  test  running  was  done  to  demonstrate  that  a  new 
and  workable  facility  has  heen  commissioned.  The  calibration  was  conducted  in  three  phases: 

(i)  Empty  tunnel  wall  pressure  distribution  measurements.. 

(ii)  Model  pressure  measurements  on  tU-in  and  15-in  chord  NACA  64A410  airfoils. 

(ill)  Balance  measurements  on  the  same  models. 

The  calibration  results  are  currently  being  analyzed.  The  mam  conclusions  reached  so  far  are: 

1.  Good  empty  tunnel  wall  pressure  distributions  for  0. 15  <  M  <  1. 10. 

2.  Diffuser  second  throat  provides  adequate  M-control, 

3.  Side  wall  boundary-layer  displacement  thickness  6*  -  0. 15"  (at  a  station  15  Inches  ahead  of  balance  centre 
lines). 

4.  The  side  wall  boundary  layer  suction  arrangement  can  provide  good  stabilizing  control  over  the  wing  flow 
adjacent  to  the  walls.  However,  the  present  porous  material  produces  an  unnecessary  large  pressure  drop, 
which  reduces  the  useful  M-Re  range  of  the  facility. 

5.  The  model  pitch  drive  and  sidewall  balance  arrangement  is  working  well. 

6.  Model  measurements  Indicate  that  the  present  wall  restraints  are  close  to  open  jet. 

7.  During  high  Reynolds  number  conditions  the  Initial  thermal  shock  during  tunnel  start-up  Is  large  enough  to  pro¬ 
duce  significant  M-variation  during  a  run. 

8.  The  open  jet  emerging  from  the  Insert  sidewalls  Into  Ine  downstream  diffuser  section  produces  severe  and  un¬ 
acceptable  vibrations  in  the  main  strut-model  support  system  used  for  driving  a  traversing  wake  probe. 

After  the  present  series  of  tests  in  the  2D  insert  arc  completed  in  January  1970,,  the  facility  will  be  dismantled  to  allow 
access  to  the  full  5-ft  x  5-ft  transonic  test  section  for  other  tests.  The  insert  will  not  be  installed  again  until  the  fall 
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of  1970  but  will  then  Include  modifications  in  its  design  aimed  at  alleviating  some  of  the  problems  encountered  during  its 
first  test  period. 


T1IEOHETICA !,  TWO- DIMENSIONAL  TllANSONIC  FLOW  STUDIES 

Computer  programs  for  subcrltieal  flow  calculation  (Sells  and  Cathcrall  methods)  have  been  adapted  to  IBM  360  system. 
Pressure  distribution  for  several  series  of  airfoils  have  been  calculated,  both  for  the  calibration  of  the  transonic  insert 
(profile  NACA  64A410)  and  for  theoretical  analysis  of  several  contemporary  airfoils.  A  study  of  accuracy  of  the  original 
Sells  method  has  been  performed,  some  improvements  have  been  introduced  and  some  will  be  added  in  the  near  future. 

The  results  of  the  "exact  in  mathematical  formulae  m”  Sells  method  have  been  compared  with  another  exact  method 
(Nieuwland's).  It  appears  that  for  many  quasi-elliptical  airfoils  there  arc  significant  differences  in  pressure  distribu¬ 
tions  (near  the  peak).  Another  accuracy  test  was  made  using  circular  arc  airfoil  and  llayasi  (Ixjckheed-Georgia)  results. 
The  Cathcrall  method  (Incompressible  flow)  gives  perfect  agreement  with  llayasi,  while  the  standard  Sells  method  gives 
results  satisfactory  for  technical  application. 

Program  for  subcrltieal  transonic  flow  with  boundary  layer  has  been  developed  (Lock-Sclls-Powell-Wllby  method)  and 
several  cases  have  boon  calculated.  For  moderate  Heynolds  number  (20  *  10*)  the  boundary-layer  effects  in  subcrltieal 
flow  decrease  the  lift  by  a  few  percent.  However,  this  method  predicts  Incorrect  pressure  distribution  In  the  region  of 
high  subcrltieal  local  Mach  numbers  because  of  simplified  compressibility  corrections. 


WIND  TUNNEL  COItltECTIONS 

Linearized  compressible  flow  analysis  applied  to  the  study  of  the  two-dimensional  perforated  wall  Interference  at  sub¬ 
sonic  speeds  (method  proposed  by  Baldwin,  Turner,  Knochtel)  has  been  reviewed  and  completed.  Calculated  values  of 
the  interference  velocity  and  its  derivatives  at  the  position  of  the  model  are  interpreted  (according  to  Hogers)  as  correc¬ 
tions  to  stream  Mach  number,  angle  of  attack,  lift  coefficient,  drag  coefficient,  and  moment  coefficient.  In  the  existing 
program  (written  in  Fortran  IV)  the  strengths  of  representing  singularities  arc  calculated  from  the  geometry  of  the 
model  and  the  measured  pressure-coefficient  distribution.  In  the  case  of  an  unknown  wall  porosity  parameter,  the  cor¬ 
rections  arc  computed  automatically  for  several  values  of  porosity  between  closed  wall  and  open  jet.  The  results  of 
this  study  are  being  applied  to  the  analysis  of  the  2D-lnsert  calibration  data. 


WIND  TI  NNEI.  TESTS  IN  T1IE  5- FOOT  X  5- FOOT  BLOWDOWN  WIND  TUNNEL  FOR  OUTSIDE  ORGANIZATIONS 
The  Hoeing  Company  (Seattle) 

in  November,  two-dimensional  transonic  testing  was  carried  out  on  the  first  of  four  (10-in  chord,  15-in  span)  Boeing 
airfoils  in  the  2D  transonic  Insert  of  the  NAE  G- ft  x  5- ft  wind  tunnel.  Measurements  included  model  balance  forees, 
wake  total  pressure  profiles,  and  wing  surface  pressure  distributions.  The  tests  were  run  at  M  =  0.615  -  0.725  at 
H  i  -  ti,  20,  and  35  y  10fl. 

In  order  to  supplement  the  poor  data  acquisition  capability  of  the  NAE  data  system,  Boeing  supplied  an  on-lino  analogue 
data  system  and  an  FM  tape  recording  unit  that  were  interfaced  with  the  wind  tunnel’s  regular  data  system.  In  addition, 
a  controller  unit  was  provided  by  Boeing  to  monitor  and  sequence  nil  functions  of  the  data  acquisition  process,  including 
the  wing  pitch- pause  program,  wing  pressure  scans,  wake  traverses,  and  data  taking  signals.  Owing  to  problems  asso¬ 
ciated  with  interfacing  and  controller  operations,  as  well  as  mechanical  breakdowns  in  the  wind  tunnel,  only  a  fraction 
oi  the  intended  program  was  completed. 

Further  testing  of  the  Boeing  airfoils  will  probably  be  continued  in  late  1070. 

SAAB  Sweden 

A  series  of  tests  involving  force  and  pressure  measurements  on  various  airfoil  configurations  have  been  carried  out  in 
the  21)  transonic  insert.  The  test  program  included  flapped  airfoils  at  low  subsonic  speed  as  well  as  "clean"  airfoils  at 
transonic  speed. 


HYDRAULICS  LABORATORY 


ST.  IAWRKNCE  SHIP  CHANNEL 

Under  the  sponsorship  of  the  Department  of  Transport,  a  study  oi  navigation,  water  levels,  tide  and  lee  problems. 
WAVE  RECORDER  DF.VEIOPMENT 

Development  of  a  staff  gauge  wave  recorder,  cither  to  lie  mounted  io  a  fixed  mast  In  shallow  water  or  floating  In  deep 
water.  Evaluation  of  other  types  of  wave  recorders:  accelerometers,  pressure  transducers,  etc, 

WAVE  INTERACTION  STUDY 

An  experimental  Investigation  to  study  the  development  of  n  wave  cascade  resulting  from  the  dynnmlc  Instability  of  a 
uniform  progressive  wave  train  of  finite  amplitude. 
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WAVE  DIRECTION  STUDY 

Field  Investigation  to  study  the  direction  of  propagation  of  wave  energy  In  a  confused  sea  as  a  function  of  wind  direction. 


IjOCK  study 

General  study  of  navigation  locks  and  locking  systems,  with  particular  emphasis  on  winter  operation. 


SEDIMENT  PROBLEMS 

To  define  and  analyze  sediment  problems  In  navigable  waters. 


VTSANUAPATNAN  HARBOUR,  INDIA 

Model  studies  to  find  remedial  measures  for  silting  and  beach  erosion. 


f/VSTRL/ME NTS  LABORATORY 


WAVE-MEASURING  DEVICES 

Modification  of  circuitry  of  absolute  type  wave  sensor  to  accommodate  a  commercially  available  form  of  voltage- 
controlled  oscillator. 


ST.  LAWRENCE  RIVER  MODEL  CONTROL  SYSTEM 

Proof  testing  of  the  weir  and  pump  control  system  completed.  Programs  for  the  control  of  the  three  weirs  and  pump 
system  by  digital  computer  being  written. 


NON-CONTACTING  PROXIMITY  GAUGE 

Preparation  of  laboratory  report. 


AIRCRAFT  BUND  LANDING  SYSTEMS 

Proposals  for  a  simple  system,  particularly  suited  for  use  in  remote  areas,  e.g. ,  Northern  Canada,  summarized  for 
Internal  distribution. 


AIRCRAFT  NOISE  PROBLEM 

Continued  acquisition  of  components  for  a  noise  monitoring  ami  recording  system  to  be  Installed  by  the  Department  of 
Transport  at  a  major  Canadian  airport. 

Collaboration  with  the  Department  of  Transport  in  regard  to:  (1)  the  formulation  of  International  standards  on  nlrport 
noise,  (2)  the  Department's  presentation  of  the  Canadian  position  nt  the  November,  1089,  ICAO  special  meeting  on  ntr- 
craft  noise. 


COMBUSTIBILITY  STUDIES 

Development  of  technique  for  the  measurement  of  the  power  dissipated  In  nr.  exploding  wire. 


WAVE-MAKING  EQUIPMENT 

Design  of  control  system  for  synchronization  (with  variable  relative  phasing)  of  two  existing  wave  makers  to  be  used 
simultaneously  in  one  hydrnulic  model. 


MECHANICAL  ARIS  TO  SURGERY 

Evaluation  of  a  prototype  semi-automatic  staple  loader  for  the  production  form  of  the  NRC/Vogelfanger  vascular 
suturing  Instrument. 
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Development  nnd  nxxcinldy  of  componentx  for  11  hvdrmillenlly  |x.we  red  lower  limit  nr  (hex  lx. 
I.iborntory  endurance  texix  on  un  cx|x-rlmcntlil  form  of  urflflHnl  vcnmix  valve. 


I'HKKDIT  (’All  STCDIKH  AND  ('All  DYNAMICS 


Cord  Intuition  of  performance  mcnxurcmcnl  .In  loaded  enrx, 
Hu-  Inlioi-ilorv, 


of  varhmx  fornix  of  ilrnfl  itcnrx  nol  provlouxly  nvnltnldr  lo 


f".,vn|!th  proving  U'xlx  on  n  new  ('I’ll  dcxltpi  of  coni  hosier  cur. 


I'  l'Ulncci'Inn  drawing  of  n  xyxtcm  for  the  mechanical  rcitmlnt  of  newsprint  completed.  Copies  sent  lo  Ihr  CPU, 

Dcncral  (IImciihhIojih  wllli  flip  Canadian  llullwnyx,  of  the  (riirk  himtlnu  prolilom  wilh  n  view  lo  Iho  formulation  of  n  xhnred 
rpxpnrcli  program.  Axxlxtnnee  lo  Hip  ('I’ll  in  Ihp  prcimriillon  unit  Inxlronipntnllon  of  n  caboose  lo  lx>  uxcrl  for  flplil 
inciiHui'Piiipnlx, 


liililnl  evil liiiii Inn  l,v  xlimilnllon,  of  n  xpponil  cnmoicrclully  iiviiIIiiIiIp  linpncl  poxlilonpr  for  (ruin  running 


N'A  Vli  1 A  ’I’ll  >.V  KVSTKM  TKSTIN.'i 


(  onl Imr. 'il  development  of  on  on-llnp  eompotcr-lKixrxl  xyxtcm  for  In— fl IkIiI  leHtlng  nnd  control  of  nnvlKntlpn  systems. 
Design  noil  ponHtniptlon  of  associated  Interface  Inxtiiinipnlnllon  nml  software. 


LOW  SPEED  AERODYNAMICS  LABORATORY 


AKIIONAI'TICAI,  WIN'D  Tt/NNKI  TKNTK  loll  Ol'TSIDK  OllCANIZA  CIONS 

Texts  wore  piii'rlcd  ool  In  flu*  <l-ll  •  li-ll  liorir.onfnl  wind  lonncl  for  dpllnvlllnml  Alrcrnfl  of  Canada,  f.lmllpd, 


N( iN-Al'ilK INAI'TICA  I.  WIND  TCNNI'I.  AND  WATKII  Tl'NNK).  TK.N’S 

I'ori'p  iiipiiMiirpiiipnlx  lor  xpvprnl  propoxed  lirldgc  xpcilonx  for  Hip  new  liiii  rnrd  inlet  proxxlng  were  poniplolpd  In  the* 
verllml  wind  liinncl  nml  n  program  of  dvnnmli  Icxlx  wax  Ix'gun, 


nn:  Niic  .'in-  cihit  v/sioi,  wind  ttnnki. 

Seven  inn  lor  construction  nml  equipment  supply  contracts  have  pxxpnllnlly  Ix-cn  completed,  und  poninilxxfnnlnK  nnd 
pnlllirnllim  of  llii'  fnclllllpx  lx  In  progress.  'Hip  I  none  I  hnx  Ix-en  operated  up  lo  full  xpprd  nml  prellmlnnry  aerodynamic 
calibrations  were  oliininoil, 


ri, nines 


Development  work  on  fluidic  veloellv  mpoxoi'm  for  oeeanograpld  r  rexenreh  hnx  xlnrled. 


nil':  AciioDVNAMies  tie  I'Iioi»i:i.i.kii-diiivi:n  v/stoi.  i.imNts  hvsikms 

A  light -weight  model  eonxlxl lug  of  n  wing  In  ll«'  xllpxli'enin  of  n  single  pro|x-ller  lx  liclng  tixtxl  to  delcrntlne  (he  (lower 
xpeelrum  of  rexponxe  lo  artificially  generated  wlml  limmd  lurlnilenep,  Strain-gauge  linlnneex  meuxiire  five  eomponenlx 
of  p  ro|  h  *  1 1 1  ’(*  force  nnd  moment,  nml  wing  UK,  drag,  nml  rolling  moment. 


WIND  TI'NNKI,  SIMI  UTION  Of'  SI'ltl'ACi:  WIND  MTHI'CTI'IIK 

l.vperlnienlni  eompni  Ixon  In  tin*  ft- ft  pilot  wind  lonncl,  of  different  melhedx  of  generating  thick  Ixamdnrv  layers 
xlmllnr  In  xlienr  nml  Mirlmlcncc  clinrnclcrlxllcx  lo  mdornl  wlmlx  In  the  lower  I  non  fl  of  the  nlmoxphere, 


Cl  IK  KXTKMNA  I.I.Y  IIMlWN  JKT  I'lAI' 


A  prolccl  Iiiim  liccn  Inlllnh'd,  liolli  In  Hn>  tt-fl  •  !<-fl  nnd  flow  visualization  wider  Innnclx,  to  InvextlKnle  the  character¬ 
istics  of  die  externally  blown  |cl  flnp,  Thin  employs  (lie  exhaust  from  a  simulated  fnn-M  engine,  which  lx  illrecled 
lhrixi|>li  the  riup  assembly  of  n  (yplciil  dixiblc-xlotlrsl,  flnpixul  nerofoll.  The  rcxultlnK  flow  Interaction  prodilrex  high 
locnl  lift  cocfl'lclenfx,  nnd  large  deflection  angles.  The  wing  nuxlelx,  which  nre  under  eonxtroellon,  nre  liolli  |wo- 
dimenxlonnl,  lull  with  n  single  Hltntilnlcd  engine  micelle. 


low 
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LOW  TEMPERATURE  LABORATORY 


TKMI’HKA  li  lt K  PIK HII.KMK  IN  IIAII, WAY  OPKIIA  THINK 

An"  Iv.l  l«-:i  I  mi, |  ttK'tilti  1  work,  onnthielod  under  (lie  misplorH  of  I  lie  AmmooIiiIi-  . . .  KmIIwmv 

,  """V1  Inelndlnu  I  ho  low  teiii|>ei„liue  performmioe  or  „li-  limko  hvnIoiiin,  nfterrooler, 

..  ,c"  development.  mi  Invesl  Iput Ion  Into  mil  MWlleh  oiiilruiu  lloiiH  under  nevere  HIiiiiiIIo  ooii.IMIomh,  oviilimtlon  of 
win, him  mil  hw I loh  hontor  hvhIoiiim.  MoIIioIh  of  e, Kiting  oonf Iiiikhih  welded  mil  lo  oliliiln  n  hiiIIiiIiIo  loiniioriilitro  for  mil 
ln.vln({  mill  mielinrlilK  mo  Udllg  Mluilloil. 


lli:i.l(’()l>Ti:H  l)K  ICINfi 

A  study  of  hellrnptcr  lolng  prnleetlnn  Involving  llio  oviiliiiitlnn  of  vhi'Iimim  nvkIoiiim  . . ml,  fluid,  mid  Helf-Hheddlng 

iiiiiIit In Im>  Mini  llio  dovolopinonl  o|  de-lelng  ooiitrol  HynloniM  liudiidliin  loo  iloloolom, 

AIIICHA  I'T  INHTHI  \\l UNTATK IN 

I'hp  Invest  Ignllnti  ol  po mm II do  iiiodpM  of  fiiiloi'o  for  idromfl  pilot  IioikIm  imdor  lelng  n iv I  snow  oondlllopH, 

M  ISC  !•:  I.IANKt  M'S  1C  I  Nil  INVKHTIC.AIIONS 

Ann  lyt  leu  I  mid  exporlnienlnl  Invest  Ignl  Ions  of  „  non-riHilino  ntiluro,  nnd  the  Invosllgntlnn  of  ooi'lnln  iihim'oIm  of  lolng 
hI niiiltii Ion  mid  inotiKu I'oinonl , 


THAWI.Klt  ICINfi 

In  oollMl»>mll,m  wllh  llio  Department  of  Transport,  mu  Invest  Ignl  Ion  of  tho  lelng  of  fishing  trawlers  nnd  other  vessels 
I'lidor  conditions  of  fronting  son  spray,  mid  of  , Methods  of  rnmlmtllng  the  prnliloiii. 


I / IW  ’IKMI'IIItATI’UI':  APPLICATIONS  IN  MKOICAL  i-NfilNKCHINC. 

An  Investigation  Into  iHHindnry-lnyor  Ihontml  flow  nioiisiironionl  In  pulHidllo  non- Newtonian  flow. 

Design  of  n  Mhuni  vnlvo  for  ooroliroMpInnl  fluid  In  hydroeephnluM,  Tho  dovloo  run  lie  ndlnsled  nfler  Implnnfnllon  without 
surglenl  Intervention. 

Development  of  t,  rnpsMlo  for  telemetry  of  Intmoriinlitl  proMMiiro.  A  Mpoolnl  fotdnro  Im  IIm  iM'rmnnonl,  rndlo-lrnnsparenl, 
in, iIhIi, ro- 1 nipcriiu'Mlde  onttlnif. 


SHIP  LABORATORY 


CAULK  PA HUNCH 

Driijf  inoMMiiroinontH  oontlnuod  on  various  configurations. 

WKATIIKIIHIIIP  SKAKKKPINO 

VnflouM  model  nntl-pltehlng  fins  nnd  Isilhs  lonlod, 

WAVk  AND  KTKKSK  ANALYSIS 

Further  work  wiim  onrrloil  iwt  In  Improvinit  existing  met  ho.  In  and  developing  now  oloolfonlo  mid  oompnlor  IooIiiiI(|iioh, 
A  nmnlH'i’  ol  Mt'ti  roeordn  wore  pnM’t'itMod, 


TANKKIf 

Now  model  prppm'cd  for  sloe  ring  assessment,  llmlio-oonlrollod  manoeuvring  ex|sTloienl,  eon, plolod  with  Iwoslerneonflg 
u rations. 

Iff:  IHIKAKINd  IN  AIICTIC 

OliMoi'ViillonK  nnd  Homo  measurements  lnkononlioi,rdCiimidliin  Icebreaker  ueeompnn.vlng  "MmihnUim"  now-being  analyzed, 
YACHT  DKHION 

llio  effort  of  various  parameters  on  the  porfornninee  of  It, rue  mil  linn  yiiehlH  nre  In-Inn  Investigated  tlioorollonlly  II 
Im  I'Xpoolod  Hint  n  limited  experimental  |iroffrnni  will  ronull. 


I 


STRUCTURES  AND  MATERIALS  LABORATORY 


i'Atici  i-:  ok  mktai.s 

«t*-«-nKth  of  welded  nuiniKfnK  Ktcel  plate;  eharnelerl*tlc*  of  atructurnl  (trade  *loel  boll*  under  fallout-  londlnit; 
1  * 1  M  ol  environment  on  lullitue  Ntrcnitth  of  II  hi  alloy*;  development  of  lout  method*  for  variable  amplitude  load* 
Mini  InvcKlIunllon  ol  effect*  of  xuch  load*;  certification  tc*t«  on  nlrrrnfl  enmirancnt*. 


IIKSPONSK  ok  KTltDCTCHKS  TO  IIKill  INTKNS1TY  NOISK 

Sluily  of  excitation  nnd  Hlruetine  rc*pon*c  merhunlnm*;  *ludy  of  panel  damplnn  ehurnc'erlatlc*  and  critical  re*pon»o 
mode*;  Invi-KllKiitlon  of  fnllKUc  dm. nine  law*;  Induttrtnl  hardware  I'ViiUinlloni  InvcKtlnatlon  of  let  cxhuuxt  notno, 


iiioii  ri ; m i>i:ha  1 1  it i-:  aiuciiakt  iiydkaki.ics 

Simlv  of  Imilt  modulo*  of  livilrmillr  fluid*,  Includlnit  emu  I* Ion*  uni!  mixture*  In  rclnllon  to  temperature  nnil  prcxxurei 
Kindle*  of  will  performance  mill  liluh  pre**urc,  lil^h  velocity  Jet  phenomena. 


UPKIIATIUNAI.  UiADS  AND  I  .IKK  OK  AIUCIIAKT  S' Tltl'CTIHKS 


liiMliuini  nliillon  of  nlrrrnfl  lor  llir  mrnm.rriiirnl  of  flight  lonil*',  nnd  acceleration*;  fatlyuc  life  monltorlnn  mvl  analyxl* 
ol  load  nml  nrrrlrrnllon  *pcctrn;  full-*ci>lc  full kxii -  M|M'rlriini texting  of  airframe*  nnil  component*;  InvcKtlKntlon  of  rough 
ground  oprrnlloiiH. 


KIlAC  TI'IIK  Ml'K'IIANK’S  AND  KIlAf "I'OlOf IY 

lnvrnllmilloii  of  orpek  lorinnllon  nml  propnunllon  In  gradient  mIi-ckm  field*;  correlative  electron  mleroHCopo  Htudle*  of 
li'ni'lnre  atirlarni;  notch  onxlllvltv  nml  frnrlure  loughne**;  rcHlnliinee  of  mulcr l» Im  (o  high  xtrnln  rate  nnd  Impact 
loading. 


IIKSKAIICII  UN  l’lU  i  l  l  i  flVi:  ('<  lA TINl IS  KOIt  ItKKHACTOHY  MKTAIJf 

Invcxilgnllor  of  muting  rompiNind*  for  prolecllon  of  refractory  meinl  xubxtrnlc*  nl  high  temperature* l  methnda  nml 
teelmli|iii>H  ol  ronlliiti  ilepoKllloiu  Hluilv  of  Interface  dlffuMbm  rate*  nml  product*)  development  of  method*  of  evaluation 
ol  plivrilrnl  pro|ici*llc*  of  coaled  lent  coupon*. 


MKCIIANICS  AND  TIIKOllY  UK  STItl'CTI  ItKS 

I’rcpnrallon  of  forimilnc  ami  cliarla  for  vllirnllona  of  an  cliiallcally  aiiapcndcd  rigid  l**lyi  atudy  of  coupled  flexurnl- 
liirMlonal  vllirallon  of  n  uniform  cantilever  Ix-ani;  devidopment  of  irnnamlHHlon  matrlceaf  extension  of  lie  am  theory  to 
Inrlinle  effcclK  of  Irannverac  curvnlurci  nlinly  of  dynamic*  of  a  cattle  under  lran*vcr*e  Impact. 


CUNSTHIT  ritl.N  AND  DKVK If IPMKNT  OK  I  t.IC.HT  IMPACT  SIMt  [ATOIt 

Dcvclopmcnl  nlCanadlnn  lacllllyi  analy*l*  of  rial  Iona  I  need*  ami  co-ordlnallon  of  ullllrnll'ini  hardware  evalunllona, 


CAUIUtATION  AND  CKItTIKICATION  UK  KOHCK  MKANI  HINt!  HKVICK8 


Kaclllllc*  available  for  llic  calibration  of  government,  unlverally,  and  Induxtrlnl  equipment  Include  deadwclttht  force 
Nlnndnnl*  up  lo  loa.oao  lb,  back-to-back  calibration  of  ncecleronietera,  and  limited  calibration  of  flud  preaaurc-type 
triin*diicer». 


COMPUSKI  i:  AND  NON-MKTA  I.I.IC  MATKIIIA1.8 

Sluille*  of  re*ln*,  <To**llnklng  I'ompound*,  and  polymerization  Inltlntorai  mnlerlnl  pro|*'rlle*  and  u*e*  of  polymer*  nml 
reinforcement*  for  i'ompo*lte*i  procedure*  for  application  and  fabrlcntloni  *truclurnl  efflelencle*. 


KINITK  I’l  1.KMKNT  MK THUDS 

Development  of  refined  finite  element*  for  plate  Ix-ndlng  nml  plane  *tre**i  development  ami  appllcallon  of  a  cylindrical 
hIicII  finite  element;  application  of  finite  element*  to  flutter  problem*;  development  of  a  Kcnernl  Irlnnitulnr  nhell 
element. 
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STRUCTURES  AND  MATERIALS  LABORATORY 
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UNSTEADY  AERODYNAMICS  LABORATORY 


STUDY  OF  VISCOUS  HYPERSONIC  FLOW  ON  A  STEADY  PLAIT. 

Surface  pressure  distributions  and  flow  field  su-veys  at  hypersonic  Mach  numbers  on  a  flat  pEitt  at  incidence,  including 
the  use  of  anemometry  and  glow-discharge  techniques 


STUDY  OF  OSCILLATORY  PHENOMENA  IN  HYPERSONIC  Flow 

Theoretical  and  experimental  investigation  of  the  effects  of  viscosity  ami  leading-edge  bluntness  on  aerodynamic  charac¬ 
teristics  of  two-dimensional  oscillating  aerofoil-  Pressure  distribution  and  flow  field  measurements;  de'ermination  of 
static  and  dynamic  stability  derivatives;  optical  studies.  Ail  experiments  performed  in  helium  at  Mach  ouml<eis  !>  and 
18.. 


PRESSURE  MEASUREMENTS  ON  OSCILLATING  WINGS  IN  SI  PERSONIC  FLOW 

Measurements  of  the  mean  surface  pressure  and  of  the  amplitude  and  phase  of  the  oscillatory  pressure  on  a  3  percent 
thick  delta  wing  at  Macii  number  1 .  *)  in  air 


HELIUM  F  YPERSONIC  WIND  TINNEL 

Two  11-inch  diameter  contoured  nozzles  (on  long  term  loan  from  the  l  S. )  for  Mach  number.-  10  and  is  now  available 
Also  available,  a  healer  to  provide  small  1 100 ’ I')  increases  in  the  stagnation  t<  mpcraturc  ol  the  flow,  and  a  vat  table 
incidence  sector  support  for  static  force,  moment,  and  pressure  measurements  on  -ting- mounted  modeLs 


DYNAMIC  EXPERIMENTS  ON  ELLIPTIC  CONES 

Measurement  of  oscillatory  characteristics  of  a  series  of  elliptic  corns  at  Mach  number  11  in  helium 


AXISYMMFTRIC  FREE  JETS 

A  method-of-characteristics  solution  of  the  flow  field  in  the  core  of  a  free  helium  let  issuing  into  a  vacuum 


SUBLIMATION  STUDIES 

Study  of  surface  recession  rates  and  of  boundary  layer  pitot  profiles  on  carbon-dioxide  flat-plate  models  in  the  helium 
wind  tunnel. 


ANALYTICAL  STUDIES  OF  UNSTEADY  FLOWS 

Analysis  of  invlscid  unsteady  flow  fields  over  elliptic  cones  m  supersonic  flows  and  over  parabolic-arc  aerofoils  in 
hypersonic  (lows. 


GAS  PHASE  REACTION  KINETICS 

A  physico-chemical  study  of  the  reactions  of  some  constituents  ol  the  upper  atmosphere,  such  as  oxygen,  mtiogen, 
and  hydrogen  atoms,  to  provide  laboratory  data  in  support  of  rocket  experiments  on  upper  atmospheric  composition, 
Proposal  available  for  a  new  method  for  the  absolute  determination  of  the  atomic  concentration  of  oxygen  and  nitrogen 
in  the  90-120  km  region. 


COSMIC  RADIATION 

Review  of  possible  hazards  for  high  altitude  flight  owing  tc  cosmic  radiation. 


PUBLICATIONS 


The  following  unclassified  reports  were  released  during  the  quarter; 


AERONAUTICAL  REPORTS 


LR-529  APPUCATION  OF  THE  HIGH  PRECISION  ThlANGULAR  PLATE-BENDING  ELEMENT  TO  PROBIEMS  WITH 

CURVED  BOUNDARIES. 

M.W.  Chernuka,  G.R.  Cowper,  G.M.  IJndberg  and  M.D.  Olson,  National  Aeronautical  Establishment, 

October  1969. 

The  application  of  the  high  precision  triangular  plate-bending  element  to  problem.-  with  curved  boundaries  is  considered. 
Appropriate  edge  conditions  for  nodal  points  on  these  boundaries  are  derived  The  error  inherent  in  representing  a  curved 
boundary  by  a  series  of  straight  segments  is  found  to  be  the  limiting  factor  on  accuracy,  while  the  effect  of  approximations 
in  the  boundary  conditions  is  minor.  To  overcome  the  first  type  of  error,  the  high  precision  element  is  modified  to  include 
one  curved  edge.  Substantial  improvements  in  accuracy  are  obtained,  as  demonstrated  in  example  calculations  for  circular 
and  elliptical  plates. 


MECHANICAL  ENGINEERING  REPORTS 


ME-232  A  TECHNIQUE  FOR  MEASURING,  RECORDING  AND  ANALYZING  TORSIONAL  OSCILI-ATTONS  OF  CRANK- 

S HAFTS  OF  RECIPROCATING  ENGINES 

r.H  liammell  and  J  E  Carvtsh,  Division  of  Mechanical  Engineering  June  1969. 

Tests  were  made  on  three  types  of  transducers,  and  a  tamable  tapai  itor  unit  that  frequency  modulates  a  carrier  wave  was 
considered  best  for  general  use  The  inlierently  good  signal-to-noise  ratio,  plus  the  ability  to  calibrate  statically  in  the 
field,  were  considered  particularly  attractive 

It  was  found  that  dynamic  calibration  of  the  transducer  should  be  carried  out  periodically,  and  a  good  Hookes  couple  is 
ideal  for  trus  purpose. 

Several  instrumentation  cheek  methods  are  outlined  that  may  be  followed  to  assure  consistency  in  the  data  produced. 


ML-4  SLEEP  DEPRIVATION  EFFECTS  ON  RESPONSES  TO  SIGNALS  OF  UNEQUAL  PROBABILITY  IN  AN  ARTI¬ 

FICIALLY  ELECTRICALLY  CHARGED  ENVIRONMENT 

L  Buck,  C  B.  Gibbs  and  R  U'onardo,  Division  of  Mechanical  Engineering,  October  1969, 

Young  male  subjects  made  repealed  runs  on  a  step-tracking  task  over  a  period  of  48  hours  during  which  they  were  deprived 
of  sleep.  Successive  runs  showed  a  reaction  time  (but  not  error  rate)  increase  that  was  inversely  related  to  signal  prob¬ 
ability,  and  a  movement  time  increase  that  was  related  to  distance  travelled.  Both  effects  were  more  evident  on  machine- 
paced  than  subject-paced  tracking.  There  was  no  evidence  that  performance  deterioration  could  be  modified  by  the  use  of 
a  patented  anti- fatigue  device. 


MP-53  DIMENSION  A  L  CHANGES  IN  SHOT  GUN  HARHELS  CAUSED  BY  THE  FIRING  OF  HARD  METAL  PELLETS. 

C.  Dayson  and  T.  Maloney,  Division  of  Mechanical  Engineering,  July  19G9. 

This  report  describes  progress  to  date  on  the  development  of  a  means  of  assessing  the  wear  and  deformation  of  shot  gun 
barrels  caused  by  the  firing  of  non-toxic  but  relatively  nard  metal  pellets.  The  results  of  an  initial  series  o'  tests  on  hard 
metal  pellets,  some  of  which  were  coated  with  non-metallle  anti-wear  coatings,  are  also  presented.  These  indicate  that 
the  enlargement  of  the  bore  of  the  barrel  at  the  nturzle  constriction  (choke)  was,  in  this  case,  primarily  due  to  the 
"hammering  out"  action  of  the  hard  pellets,  with  the  wear  of  the  bore  lieing  a  secondary  etfect.  The  pellet  coatings  tested 
proved  to  be  Ineffective  as  a  means  of  preventing  either  form  of  damage, 


MS-119  AUTOMOBILE  DEFROSTER  PERFORMANCE  DURING  A  CANADIAN  WINTER. 

G.F.W.  McCaffrey,  National  Aeronautical  Establishment,  June  1968. 

Considerable  dissatisfaction  has  been  expressed,  by  drivers,  regarding  the  performance  of  defroster  systems  in  Canadian 
winter  driving  conditions.  In  an  attempt  to  obtain  some  data  on  this  problem  a  group  of  six  typical  automobiles  was 
subjected  to  a  standardized  test  of  defroster  performance  and,  subsequently,  monitored  during  normal  driving  from 
November  1967  to  March  1968. 
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MECHANICAL  KNC.LNEEIHNG  REPORTS  (Cunt'd) 


MS  US  A  Cr-TI-SI  DIFFUSION  COAI1NG  KOH  THE  NIOBIUM  ALLOY  BOO:  AN  INVESTIGATION  OF  SOME  COATING 

PAHA  METERS. 

J.M  Trenouth  and  M.S.  GrzcdzieLski,  National  Aeronautical  Establishment.  Dctt  ml*  r  1909 


Chromium-titanium  and  chromium-titanlum-silicoo  diffusion  routings  Here  formed  on  smeimens  of  the  niobium  alloy  B66 
by  vacuum  pack  cementation.  Some  of  the  coating  parameters  for  formation  of  the  chromiuin-titanium  diffusion  coaling 
were  investigated  with  respect  to  the  amount  and  structure  of  the  coating  formed.  The  rate  i fei  mation  of  this  coating 
was  found  to  be  diftusioa-dependent.  Consequently  the  rale  was  strongly  dependent  on  the  processing  parameters  of  time 
and  temperature,  and  essentially  independent  of  vapour-phase  conditions  such  as  activator  concentration  and  pressure 
(evacuation  rate).  The  rate  of  formation,  and  the  structure  and  composition  of  tlic  icsulting  coating  were  a  function  of  the 
pack  alloy  and  hence  of  the  coating  atmosphere  composition.  For  coatings  formed  in  a  60/40  w,o  Cr/Tt  alloy  pack,  the 


overall  formation  rate  was  calculated  to  be  KCl  .,  -  2. 14  /  10°  exp  (  •-  ^ 

two  zones  of  the  coating  (overlay  and  diffusion  zone)  were  determined  to  be  K0,,lUj 
1.5s  x  105  exp  (-  ) 


.I'l 
p  h 

1.07 


Thi*  lot  mation  rates  for  the 
42  >  102' 


10  exp 


/  42  ■  ion  .1,1 

V  HT  /  M/h 


and  K 


Diffupcn  ’on* 


gh  .  The  amount  of  diffusion  alloying  of  silicon  with  previously 


formed  chromium-titanium  coatings  was  indepewient  of  the  chromium-titanium  coaling  conditions  and  thickness  As  a 
result  of  this  constant  silicklc  layer  thickness  on  different  thicknesses  of  the  chromium-titanium  coating,  van  ing  structures 
of  the  final  chromlum-tltanium-sHtcon  coating  were  obtained. 


MS-127  THE  INFLl/ENCE  OF  SODIUM  SULPHATE  AND  SODll  M  CHLORIDE  ON  THE  OXIDATION  RESISTANCE  OF 

A  Cr- 11-Si  COATING  ON  ALLOY  B66. 

E.P.  Whelan  and  J.  M.  Trenouth,  National  Aeronautical  Establishment,  December  1969., 

An  investigation  has  been  made  of  the  influence  of  Na;S04  and  NaCl  on  the  oxidation  resistance  of  Cr-Ti-Si  coatings  de¬ 
posited  by  pack  cementation  on  the  niobium  alloy  B66. 

Coated  specimens  were  contaminated  by  Na;i>04  salt  concentrations  of  approximately  1  Mg/c>n\  and  oxidised  at  11  litre, 
1200*0,  and  1300°C.  Oxidised  contamlnalcd  specimens  exhibited  initial  mass  gains  considerably  greater  than  those  of 
identical  imeontaminated  specimens.  The  contaminated  surface  lieeame  brown,  glassy,  aud  bubbly  within  io  minutes  of 
being  exposed  to  the  oxidising  environment. 

Microstruetures  of  oxidised  contaminated  specimens  indicated  that  the  protective  silicidc  coating  zones  oxidised  at  a  more 
rapid  rate  than  did  those  of  uncontaminated  specimens,  and  that  an  array  oi  particles,  possibly  titanium  oxides  or  nitrides, 
formed  in  the  diffusion  zone  in  less  than  i00  minutes  at  1300*C.  Electron  probe  microanalysis  did  not  detect  any  sodium  or 
sulphur  In  the  oxide  or  substrate  after  300  minutes  exposure  at  1300"C.  Oxidation  prior  to  Na2S04  contamination  did  not 
increase  the  subsequent  coi  rosion  resistance  of  the  coatings 

A  preliminary  investigation  indicated  tnat  NaCl  also  influenced  the  oxidation  resistance  of  coated  specimens  in  a  deleteri¬ 
ous  manner. 


MISCELLANEOUS  PAPERS 


BUCK,  L  How  Does  Drinking  Impair  Driving’  (Comment  l'alcool  sffaibiit  las  refitxes  au  volant).  Royal  Canadian  Mounted 
Police  Gazette,  31(10),  18-22,  October  I960  (English  and  French  editions) 

CHAN,  Y.  Y.  Parametric  Approximation  for  Incompressible  Laminar  Boundary  layers  with  Suction  ant1  Injection.  CASI  Trans¬ 
actions  Vo)  2,  No.  2,  September  1369,. 

CHAN,  Y.y.  An  Experimental  Study  cl  a  Yawed  Circular  Cone  In  Hypersonic  Flows  AIAA  Journal,  Vol  7,  No.  10,  October  I960 

COWPER,- G.  R. ,  KOSKO,  E,,  LJNDBERG,  G.  M.  and  OLSON,  M.D  Static  and  Dynamic  Applications  of  a  High-P.’t. vision  Plate 
Bending  Element.  AIAA  Journal,  Vol,  7,  No.  10,  October  1969,  pp,  1957-1965 

DEBLOIS,  J.  D  Changing  Time  -  Cinematography  at  Speeds  Other  Than  Norma!..  To  be  published  in  the  Proceedings  of  the  1969 
Symposium  of  the  Canadian  Science  Film  Association, 

FIS2DON,  W. ,  IYENGAR,  S.  and  ORL1K-HUC KE MANN,  K  J.  On  Unsteady  Effects  of  Viscosity  on  the  Hypersonic  Flow  Past 
a  Slowly  Oscillating  Wedge  with  Attached  Shock  Wave.  National  Research  Council  of  Canada,  NAE  Lab,  Tech.  Report 
LTft-UA-11,;  September  1969, 

GARDNER,  1  Colorimetnc  Determination  of  Ethylene  Glycol  Mononiethyl  Ether  In  Aviation  Turbine  Fuel,  Journal  of  the 
Institute  of  Petroleum,  Vol,  55,  No.  546,  November  1969, 

LEE,  B..M.K  A  Modified  Shock  Layer  Theory  tor  Hypersonic  Internal  Flows,  CASI  Transnit.'ons  Vol,  2,.  No.  2,  September 
1969,  pp.  67-74. 

LEE,  B.  H.K  The  Pi  eduction  of  Strong  Shocks  in  a  Multi-Stage  Gaseous  Detonation  Driven  Shock  tube.  CASI  Transactions 
Vol,  2,  No,  2.  September  1969,  pp.  75-80. 


MEYER,  H  F,  and  GRAHAM,  G  D.  A  Non-Axisymmotrlc  Hypersonic  Blast  Wave  Analogy.  CASI  Transactions, 
September  1969. 
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MISCELLANEOUS  PAPERS  ICat’dl 


MUFTI,  1.  H.  CHOW,  C.K.  and  STOCK,  F.  T.  Solution  of  Ill-Conditioned  Linear  Two- Point  Boundary  Value  ProMeaa  fay  the 
Riccali  T ranaformatlon .  SIAM  Review,  Vol.  11,  No.  4.  October  1969. 

ORUK-RLCKE  V  *.KN,  K  J.  Simple  Formulas  for  Unsteady  Pressure  on  Slender  Wedges  and  Cones  In  Hypersonic  Flow, 

J.  Spacecraft  and  Rockets,  Vol.  6,  No.  10,  October  1969,  pp.  1209-1211. 

PEAKE,  D.  J.  Three-Dimensional  Flow  Separations  on  Upswept  Rear  Fuselages,  CASJ  Journal,  December  1999. 

ROMERO- SIERRA,  C.,  TANNER,  J.A.  and  VILLA,  F.  EMC  Changes  in  the  Limb  Muscles  of  Chickens  Atbjected  to  Micro- 
wave  Radiation.  National  Research  Council  of  Canada,  DME  Lab.  Tech.  Report  LTR-CS-10.  November  1969. 

HOMERO-SIEKKA,  C. ,  HUSBAND,  C.  and  TANNER,  J.  A,  Effects  of  Microwave  KadUtion  on  Parakeets  in  Flight.  National 
Research  Council  of  Canada,  DME  Lab.  Tech.  Report  LTR-CS-1S,  November  1969. 

TANNER,  J.  A. ,  ROMERO- SIERRA,  C,  and  DAVIE,  S.  J..  A  Preliminary  Investigation  of  the  Effect  o i  Pulsed  Microwaves  an 
the  Feeding  Behaviour  of  Chickens,  National  Research  Council  of  Canada,  DME  Lab.  Tech.  Report  LTB-C8-20, 
December  1969. 

VAN  DAM,  W. ,  %’ANNER,  J.A  and  ROMERO-S1ERRA,  C.  A  Preliminary  Investigation  of  Piezoelectric  Effects  ill  Bird 
Feathers  National  Research  Council  rf  Canada,  DME  Lab.  Tecb.  Report  LTR-CS-19,  December  1969. 

VILLA,  F. ,  ROMERO-S1E11RA.  C.  and  T  >NNER,  J.A.  Techniques  for  Extracting  Bioelectric  Signals  from  Birds  Subjected 
to  Microwave  Radiation.  National  Research  Council  of  Canada,  DME  Lab.  Tech.  Report  LTR-CS-15,  November  1909. 

WARD  LAW,  R  L.  and  PONDER,  C.  A.  An  Example  cf  the  Use  of  Wind  Tunnels  for  Investigating  the  Aerodynamic  Stability  of 
Bridges.  Proceedings  Canadian  Good  Roads  Association  Convention,  September-October  1969. 

WUODSIDE,  C.  M  Pellet  Steelmaking  Modelling  and  Control  -  Interim  Report  No.  1,  National  Research  Council  of  Canada, 
DME  Lab,  Tech.  Report  LTH-CS-14,  December  1969. 


UNPUBLISHED  PAPERS  AND  LECTURES 


BAKER.  H  and  GRAEFE  PW.  U.  Ocean  Wave  Profiles  and  Spectra  Measurements  Using  an  Airborne  Magnetometer. 

Presented  at  Annual  Meeting  ot  the  Canadian  Committee  on  Oceanography  and  Second  International  Oceanographic  Symposia, 
Empress  Hotel,  Victoria,  B  C.,  17-21  November  1969. 

COWPER,  G.  It.  Finite  Element  Methods  for  the  Stress  Analyst.  Lecture  presented  to  the  Society  for  Experimental  Stress 
Analysis,  in  Ottawa,  Ontario,  2  December  1969, 

DAYSON,  C.  A  New  Type  of  Stepped  Thrust  Bearing.  Presented  at  Graduate  Student  Seminar,  University  of  Alberta,  Edmonton, 
Alberta,  8  September  1969. 

DEBLOIS,  J.  D.  The  Motion  Picture  Camera  in  Scientific  Investigation.  Presented  at  the  Second  Annual  Conference  on  Bio- 
Medical  Photography  at  Rochester  Institute  of  Technology,  Rochester,  New  York,  25  October  1969, 

DEBLOIS,  J.D.  High  Speed  Camera  Techniques.  Presented  at  the  Canadian  Society  of  Cinematographers,  Ottawa  Meeting, 

29  November  1969. 

GODBY,  E.A.,  HOOD,  P.  J.*  and  BOWER,  M.*  Aeromagnetic  Reconnaissance  of  Baffin  Bay,  the  Labrador  Sea  and  the  North 
Atlantic  Ocean.  Presented  at  Annual  Meeting  of  the  Canadian  Committee  on  Ocee:. -graph}-  and  Second  International 
Oceanographic  Symposia,  Empress  Hotel,  Victoria,  B,C.,:  17-21  November  1961. 

INCE,  S.  Winter  Navigation  Problems  on  the  St.  Lawrence  River.  Seminar  at  the  Ur. versify  of  Toronto,  Toronto,  Ontario, 
November  1969, 

KASVAND,  T.  Computing  Machines  and  Intelligence.  Lecture  (in  Estonian)  presented  to  the  Academic  People  (Estonians)  and 
guests  from  Latvian  and  Lithuanian  Societies  for  the  50th  Anniversary  of  the  Estonian  National  University  of  Tartu, 

1  December  1969. 

KASVAND,  T.  Computers.  Lecture  (in  Estonian)  presented  to  the  Estonian  Society  In  Ottawa,  Ontario,  1  November  1969. 

LINDBERG,  G  M.  and  OLSON,  M,D.  Development  of  High  Precision  Triangular  Plate  Bending  anJ  Shallow  Shell  Finite 
Elements.  Seminar  given  at  University  of  Waterloo,  Waterloo,  Ontario,  25  November  1969. 

LINDBERG,  G.M.  and  OLSON,  M.D,  Applications  of  Finite  Element  Techniques.  Lecture  presented  to  the  Society  for  Experi¬ 
mental  Stress  Analysis,  In  Toronto,  Ontario,  25  November  1969, 

MILNER,  M.  Studies  of  Human  locomotion  and  Programmed  Stimulation  of  Skeletal  Muscle.  Lecture  presented  to  Prosthetics 
and  Orthotics  Research  and  Development  Unit,  Manitoba  Rehabilitation  Hospital;  Research  Group,  Shriner's  Hospital  for 
Crippled  Children,  Medical  Staff,  University  of  Manitoba,  Winnipeg,  Manitoba,  7  November  1969. 

MILNER,  M.  A  Bio-Engineer  Looks  at  Huma’i  locomotion.  Lecture  presented  to  Emory  University  Regional  Rehabilitation, 
Research  and  Training  Center,  Atlanta,  Georgia,  U.S,  A  ,  4  December  1969. 


D.  E.  M.  R  Geological  Survey  of  Canada. 
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WWlMltD  PAPEB8  AMD  UOUtB  tCcct'dl 

MUB,  M.  BrniiHlw  In  TWn  Engineering  Engineering  94.57* ,  a  grmhiatr  cone  of  etndy  presented  to  Ike  Engineering 
Dapnrtnmnt  at  Carktoa  University.  Ottawa,  Ontario.  tmiati  r-Pwcabtr  1M. 

PC  LUG.  J.  Pttfcatoa  «f  **TI  hXo  Miobhun  angle  Crystals.  TMS-AP1E  Fall  Meetkg,  PhitadrlpMs.  Pa..  U.S.A.. 

Onakar  1M.  Paper  preesMedby  Dr.  E.P.  Wheta  in  the  absence  of  Or.  J.  Pelleg. 

IIEVBUDM,  H.T.  A  gwfct  and  Recent  DwriyH  It  Creeh  Poahtoa  Indicator  -  Soft  Landtag  Catena.  PicaeMed  Is  the 
Caaadtaa  AercaaaUca  aad  S^aao*  InaUtatc,  Ottawa  Drench.  Ottawa,  Oatario,  It  December  IMS. 

TAKNET,  Jf.W.  Bane  Pmhfaatt  lawhhi  the  t?ace  Ice  jJatinaadtd  Tartateat  Jets  ia  Fluidic  Sensors.  Talk  given  to  the 
staff  of  the  SACLANT  ASH’  Beaearch  C eater.  La  SperJi,  Italy,  22  September  IMS. 

WHELAN,  E.P.  ami  DAINTY,  R.  V.  The  Cahlatiaa  Resistance  of  Commercial  AhimtaMe  Coalite  oa  the  SbperaUoys  1NCO  T13C 
aral  01 1M.  Paper  presented  by  Or.  E.P.  WheJwi  at  TliS-ADIE  Fall  Meetfc*,  Philadelphia.  Pa..  I.'.  S.  A,  16  October  1969. 


AERONAUTICAL  LIBRARY 

Statistical  Summary  of  Library  Operations  for  the  Quarter 
_ October  -  December  1969 _ 


Documents  accessioned  (including  duplicates) .  3, 728 

Documents  accessioned  (Drat  copier  only) . , .  3, 329 

Car  A  added  to  the  catalogue .  13, 021 

Boohs  received . 237 

Bound  periodicals  received .  167 

Loans  to  NRC  staff  (including  Periodical  circulation  and  Xerox  and  Microfiche 

copies  supplied  in  lieu  of  loans! .  7, 765 

Loans  and  distribution  to  outsiders .  2, 339 

Total  circulation .  Id,  104 

Statistical  Summary  for  1968  and  1969 

1968  1969 

Documents  accessioned  (Including  duplicates) .  . .  15, 268  16, 886 

Documents  accessioned  (first  copies  only) . 11,461  12,495 

Cards  added  to  the  catalogue . . .  62, 853  57, 800 

Books  received .  1,557  1,587 

Bound  periodicals  received . 337  648 

loans  to  NEC  staff  (including  Periodical  circulation  and  Xerox  and  Microfiche 

copies  supplied  In  lieu  of  loans) . .  30, 505  30, 042 

Loans  and  distribution  to  outsiders . .  7, 648  9,138 

Total  circulation .  38, 153  39, 180 


NOTEi 


These  summaries  include  statistics  for  the  Uplands  Branch  of  the  Aeronautical  Library. 


93 


PROPRIETARY  PROJECTS  DURING  1969 


Put  of  the  work  of  the  two  Divisions  covers  proprietary  project!, 
and,  for  thle  reaeon,  has  not  been  reported  in  these  Bulletins. 

Following  is  a  l‘st  of  Industrial  Organizations,  Government  Depart¬ 
ments  and  Universities  for  whom  work  was  done  daring  IMS. 


INDUSTRIAL  ORGANIZATIONS 


Aero  Photo  Inc.,  Quebec,  P.Q. 

Aerospace  Exhibits  Limited,  Montreal,  P.Q. 

Air  Canada,  Montreal,  P.Q. 

Aircraft  Industries  of  Canada  Ltd. ,  St. Jean,  P.Q. 

Allatt  Auto  Supply,  Toronto,  Out. 

Automatic  Electric  Co.,  BrockviUe,  Out. 

Avian  Aircraft  Ltd.,  Georgetown,  Out. 

Babcock -Wilcox  Canada  Ltd. ,  Galt,  Ont. 

Back-Simpson  Ltd.,  London.  Ont. 

Barnes  Engineering  Co.,  Stamford,  Conn.,  U.S.A. 
Battelle  Memor.  '1  Institute,  Columbus,  Ohio 
Bell  Helicopter  Co.,  Fort  Worth,  Texas,  U.S.A. 

Bobtcx  Corporation  Ltd.,  Montreal,  P.Q. 

Boeing  Company,  Seattle,  Washington,  U.S.A. 

Boeing  of  Canada  Ltd. ,  Arnprior,  Ont. 

Bourns  (Canada)  ltd.,  Toronto,  Ont. 

B.P.  Canada  Ltd.,  Oakville,  Ont. 

B.P.  Canada  Ltd.,  Montreal,  P.Q. 

Bristol  Aerospace  Industries,  Montreal,  P.Q. 

Bristol  Aerospace  Ltd,,  Winnipeg,  Man. 

Bristow  Instruments  Company,  Edmonton,  Alta. 

CAE  Industries  Ltd.,  St. Laurent,  P.Q. 

CAMAT  Transportation  Consultants  Inc. ,  Montreal,  P.Q.. 
Canadair  Ltd.,  Montreal,  P.Q. 

Canadian  Aero  Services  Limited,  Ottawa,  Ont. 

Canadian  Amateur  Ski  Association,  Montreal,  P.Q.- 
Canadian  Aniline  &  Extract  Co.,  Hamilton,  Ont, 

Canadian  Aviation  Electronics,  Montreal,  P.Q. 

Canadian  Broadcasting  Co. ,  Ottawa,  Ont. 

Canadian  Car  Pacific,  Vancouver,  B.C. 
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Canadian  Gas  Association,  Toronto,  Ont. 

Canadian  General  Electric  Co.,  Ltd.,  St. Andrew*  East,  P.Q. 
Canadian  General  Electric,  Port  Hawkesbury,  N.S. 

Canadian  Industries  limited,  McMasterville,  P.Q. 

Canadian  Inge raol -Rand  Co. .  Montreal,  P.Q. 

Canadian  International  Paper  Company,  Gatineau.  P.Q. 

Canadian  Liquid  Air,  Montreal,  P.Q. 

Canadian  National  Railway  a.  Montreal,  P.Q. 

Canadian  Pacific  Airlines,  Vancouver,  B.C. 

Canadian  Pacific  Railway  Co.,  Montreal,  P.Q. 

Canadian  Petrofina,  Po  In  te-Aux -Trembles,  P.Q. 

Canadian  Vickers  limited,  Montreal,  P.Q. 

Canadian  Westinghouse  Co. ,  lid. ,  Hamilton,  Ont. 

Canadian  Yacht  Consortium,  Halifax,  N.S. 

Canive  lid.,  Hagersville,  Chit. 

Capital  Air  Surveyr.  Limited,  Ottawa,  Ont. 

Centennial  Auto  Parts,  Ottawa,  Ont. 

Central  Design  and  Drafting  Ltd.,  Montreal,  P.Q. 

Century  Air  A  Ground  Services,  Montreal,  P.Q. 

Cessna  Aircraft,  Wichita,  Kansas 
Champlain  Power  Products  Ltd. ,  Toronto,  Ont. 

Charbonneuu,  Mr.  G.,  L'Origlnal,  Ont. 

Cities  Service  Oil  Co. ,  Cranbury,  N.J. 

City  of  Ottawa,  Ottawa,  Ont. 

Coast-Eldridge  Professional  Services  Division,  Vancouver,  B.C. 
La  Compagnle  Photo-Air  Laurentldes,  Quebec,  P.Q. 

Computing  Devices  of  Canada  Ltd. ,  Ottawa,  Ont. 

Consolidated  Bathurst  Paper  Co,,  Grandmere,  P.Q. 

Cook  Electric  Co. ,  Chicago,  Ill.,  U.S.A. 

Cornell  Aeronautical  Laboratory  Inc,,  Buffalo,  N.Y.,  U.S.A. 

Tj.  Corrigan  U  Co.,  Hamilton,  Ont. 

G.F.  Crate  Ltd,,  Toronto,  Ont. 

Crawley  Films  Limited,  Ottawa,  Ont. 


DeHavilland  Aircraft  of  Canada  Ltd. ,  Downsvlew,  Ont, 
DeHavilland  Aircraft  of  Canada  Ltd, ,  Longueuil,  P.Q. 
Dowty  Equipment  of  Canada  Ltd,,  Ajax  Ont. 

Dufresne  Piling  Company,  Ottawa,  Ont. 

DuPont  of  Canada  Ltd.,  Montreal,  P.Q,. 

Dustbane  Ltd.,  Ottawa,  Ont. 


Eastern  Air  Lines  Inc.,  New  York,  N.Y.,  U.8.A. 

Eastman  Kodak  Co.,  Rochester,  N.Y. 

EMI  Electronics  Ltd. ,  Dartmouth,  N.8. 

Lnjay  Additives  Laboratory,  Linden,  N. J. 

Esso  Research  and  Engineering,  Linden,  N.J. 

Fathom  Oceanology  limited.  Port  Credit,  Out. 

Ferranti -Packard  Ltd.,  Toronto,  Ont. 

Field  Aviation  Ltd. ,  Malton,  Ont. 

Fleet  Manufacturing  Ltd.,  Fort  Erie,  Ont. 

Furnace  Engineering,  Montreal,  P.Q. 

Garrett Mfg.  Ltd.,  Rexdale,  Ont. 

General  Dynamics,  Convalr  Division,  San  Diego,  Calif.,  U.S.A. 

General  Electric  Co.,  Cincinnati,  Ohio 
General  Motors  of  Canada,  Ottawa,  Ont. 

German  and  Milne,  Naval  Architects  and  Marine  Surveyors,  Montreal,  P.Q. 
Golden  Eagle  Can.  Ltd.,  Montreal,  P.Q. 

Grace,  W.R.  Hatco  Div.,  Fords,  N.J. 

Great  Can.  Oil  Sands,  Fort  McMurray,  Alta. 

Gulf  Oil  Canada  Ltd. ,  Sheridan  Park,  Ont. 

Hamel,  Malouin  A  Assoc.  Druramondville,  P.Q. 

Hauts-Monts  Inc.,  Quebec,  P.Q. 

The  Hospital  for  Sick  Children,  Toronto,  Ont. 

Howe  India  Private  Limited,  New  Delhi,  India 
Hugh  Carsons,  Ottawa,  Chit. 

Huntec,  Division  of  Kentlng  Exploration  Services  Ltd. ,  Toronto,  Ont. 
Hydradrtve  Engines  Ltd.,  Toronto,  Ont. 

Hydro-Electric  Power  Commission  of  Ontario,  Toronto,  Ont. 

Ian  Martin  Associates,  Ottawa,  Ont. 

Imperial  Oil  Enterprises  Ltd.,  Charlottetown,  P.E.I., 

Imperial  Oil  Enterprises  Ltd.,  Darmouth,  N.S.: 

Imp<.  rial  Oil  Enterprises  Ltd, ,  Lewis|»rte,  Nfld. 

Imperial  Oil  Enterprises  Ltd. ,  Montreal,  P.Q. 

Imperial  Oil  Enterprises  Ltd.,  Sarnia,  Ont. 

Imperial  Oil  Limited,  Toronto,  Ont. 

Industrial  Metal  Co.  of  Canada,  Toronto,  Ont. 

International  Cellulose  Research  Ltd.,  Hawkesbury,  Ont. 

Internationa!  Paper  Co.,  Gatineau,  P.Q., 

Iron  Ore  Co.  of  Canada,  Labrador  City,  Nfld. 
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Johnson,  Matthey  &  Mallory  Ltd.,  Toronto,  Ont. 

Renting  Aviation  Limited,  Malton,  Ont. 

Knox,  Mr.  James,  Ottawa,  Ont. 

LaurentianAlr  Services  Ltd.,  Ottawa,  Ont. 

Leigh  Instruments  Ltd.,  Carleton  Place,  Ont. 

Leigh  Systems  Inc. .  Syracuse,  N.Y.,  U.S.A. 

Lockheed  Aircraft  Corp. ,  Burbank,  California.  U.S.A. 

Lockheed  Missiles  &  Space  Co. ,  Sunnyvale,  Calif. 

Lockwood  Survey  Corp.  Ltd.,  Toronto,  Ont. 

Lockwood  Su rvey  Corp .  Ltd.,  Vancouver,  B.C. 

Lubri/.ol  Corp. ,  Cleveland,  Ohio 

Machalr  Survey  Limited,  Calgary,  Alta. 

MacMlllnn  Blocricll  and  Powell  nivor  Co.,  Vancouver,  B.C. 
Mechron  Engineering  Products,  Ottawa,  Ont. 

Merck,  Sharpe  &  Dnhrne  Can.  Ltd.,  Pointe  Claire,  P.Q. 

Miltronics  Ltd.,  Peterborough,  Ont. 

Mobil  Res.  *<  Development  Co.,  Paulsboro,  N.J. 

Montreal  Transportation  Commission,  Montreal,  P.Q. 

Mulrhead  Instruments  Ltd.,  Stratford.  Ont. 

National  Aeronautics  and  Space  Administration  (USA) 

Noranda  Copper  Refining  Industries  Ltd. ,  Montreal,  P.Q. 

North  American  Rockwell  Corp. ,  Tulsa  Division,  Tulsa.  Oklahoma 
Northern  Electric  Co.  Ltd.,  Belleville,  Ont. 

Orenda  Ltd.,  Malton,  Ont. 

Ottawa  Civic  Hospital,  Ottawa,  Ont. 

Ottawa  Transportation  Commission,  Ottawa,  Ont. 

Philips  Electronics  Industries  Ltd.,  Toronto,  Oit. 

Preci-Tools  Limited,  Montreal,  P.Q. 

Public  Fuel  Transmission  Systems,  Ottawa,  Ont. 

Quebec  Cartier  Mining  Co.,  Port  Cartier,  P.Q. 

Quebec  Hardwoods  Inc.,  Hull,  P.Q. 

Quebec  North  Shore  fc  Labrador  Railway  Co.,  Sept  lies,  P.Q. 

Rails  Co.,  New  Jersey 

Range  Aerial  Survey  Limited,  Calgary.  Alta. 

Research  Engineering  Associates,  Toronto,  Ont. 

Rio  Algom  Mines  Ltd. ,  Toronto,  Ont. 
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Rolls -Hoyce  (Derby)  Ltd.,  Derby,  England 
Rolls-Royce  (Montreal)  Ltd.,  Montreal,  P.Q. 

Ron  Engineering  &  Construction  Ltd.,  Ottawa,  Ont. 

SAAR  Aktlebolag,  Sweden 

Saguenay  Terminals  Limited,  Port  Alfred,  P.Q. 

Sandia  Corporation,  Albuquerque,  New  Mexico 
Saunders  Aircraft  Corp.  Ltd.,  Dorval,  P.Q. 

Schofield  &  de  Vries,  Breslau,  Ont. 

Separator  Engineering  Ltd.,  Montreal,  P.Q. 

Shell  Canada  Ltd.,  Montreal,  P.Q. 

Shell  Canada  Ltd.,  St.  Boniface,  Manitoba 
Shell  Canada  Ltd.,  Toronto,  Ont. 

Societe  Nucletron  Inc. ,  St.  Leonard,  P.Q. 

Space  Research  Institute.  Montreal,  P.Q. 

Spar  Aerospace  Company,  Toronto,  Ont. 

Spartan  Air  Services  Ltd.,  Ottawa,  Ont. 

Sperry  Gyroscope  Company,  Montreal,  P.Q 
The  Steel  Company  of  Canada  Ltd. ,  Hamilton,  Ont. 

Sud  Aviation,  Marignanc,  France 
Sun  Oil  Canada  Ltd.,  Toronto,  Ont. 

Survair  Limited,  Ottawa,  Ont. 

Sutton  Perry  Ltd.,  Ottawa,  Ont. 

Swan  Wooster  -  CAT),  Vancouver,  B.C. 

Terry  Industries  oimited,  Dorval,  P.Q. 

Texaco  Canada  Ltd.,  Montreal,  P.Q. 

Texaco  Inc.,  Beacon,  N.Y. 

Toronto  Coppersmithing  Co.  Ltd.,  Toronto,  Ont. 

Toronto  General  Hospital,  Toronto,  Ont. 

Trans  Canada  Pipelines,  Toronto,  Ont. 

Tretolitc  Division  of  Pctrolitc  Corp. ,  Clarkson,  Ont. 
Turbo  Chemical  Co.,  Edmonton,  Al*a. 

Underwriters  Adjustment  Bureau,  Hull,  P.Q. 
Underwriters’  Laboratories  of  Canada,  Scarborough,  Ont. 
Union  Carbide  Ltd.,  Don  Mills,  Ont, 

United  Aircraft  of  Canada  Ltd. ,  Iongucuil,  P.Q. 

Vale  Enterprises,  Montreal,  P.Q. 

Vapor  Heating  Limited,  Montreal,  P.Q. 

Vnrltin  Associates  of  Canada  Ltd.,  Georgetown,  Ont. 
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Wabash  Mine*.  Point  Noir,  P.Q. 

Western  Caissons  limited,  Laval,  P.Q. 

Weatinshouae  Electric,  Philadelphia,  Pa. 

Worthington  (Canada)  lid.,  Brantford,  Onl. 

GOVERNMENT  DEPARTMENTS 

Atomic  Energy  Control  Board,  Ottawa,  Ont. 

Atomic  Energy  of  Canada  Limited,  Ottawa,  Ont. 

Atomic  Energy  of  Canada  Ltd. ,  Winnipeg,  Manitoba 

The  Bank  of  Canada,  Ottawa.  Ont. 

British  Columbia  Research  Council,  Vancouver,  R.C. 

Canada  Centre  for  Inland  Waters,  Burlington,  Cut. 

Canadian  Broadcasting  Corporation,  Ottawa,  Ont. 

Canadian  Government  Specifications  Board,  Ottawa,  Ont. 

Consumer  fc  Corporate  Affairs,  Dept,  of  Trade  fc  Practices,  Toronto,  Ont. 

Defence  Research  Board,  GARDE,  Ottawa,  Ont. 

Defence  Research  Establishment,  Halifax,  N.S. 

Defence  Research  Establishment,  Valcartier,  P.Q. 

Department  of  Agriculture,  Ottawa,  Ont. 

Department  of  Defence  Production,  Ottawa,  Ont. 

Department  of  Energy,  Mines  and  Resources,  Ottawa,  Ont. 

Department  of  Fisheries  fc  Forestry,  Ottawa,  Ont,, 

Department  of  Forestry  A  Rural  Development,  Ottawa,  Ont. 

Department  of  Indian  Affairs,  Northern  Development,  Ottawa,  Ont. 
Department  of  Industry,  Trade  A  Commerce,  Ottawa,  Ont. 

Department  of  Justice,  Ottawa,  Ont. 

Department  of  National  Defence,  Ottawa,  0-_t. 

Department  of  National  Health  U  Welfare,  Ottawa,  Ont. 

Department  of  Public  Works,  Ottawa,  Ont. 

Department  of  Transport,  Meteorological  Branch,  Toronto,  Ont. 
Department  of  Tiansport,  Ottawa,  Ont. 

Department  of  Transport,  Prescott,  Ont., 

Engineering  Research  Service,  Central  Experimental  Farm,  Ottawa,  Ont., 

Forest  Products  Laboratory,  Ottawa,  Ont. 

Hydro-Quebec,  Montreal,  P.Q. 
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Ministry  of  Technology.  Aeroplane  It  Armament  ExperlmecUl  Establishment, 

Bos  combe  Down,  Salisbury,  England 

Museum  of  Science  and  Technology,  Ottawa,  Out. 

National  Severe  Storms  Laboratory,  Norman,  Okla.,  U.S.A. 

Ontario  Department  of  Lands  &  Forests,  Ottawa.  Ont. 

Ontario  Department  of  Lands  b  Forests,  Toronto,  Ont. 

Ontario  Provincial  Government,  Gasoline  Tax  Branch,  Toronto,  Ont. 

Ontario  Provincial  Police,  Perth,  Ont. 

Research  Council  of  Alberta,  Edmonton,  .Vita. 

Royal  Canadian  Mounted  Police,  Ottawa,  Ont. 

Supreme  Court  of  Ontario,  Toronto,  Ont. 

United  States  Air  Force,  Wright -Patterson  Flight  Dynamics  Group  Dayton,  Ohio,  U.S.A. 
United  States  Army,  Paulsboro,  N.J. 

U.S.  Department  of  the  Navy,  U.S.  Naval  Air  Station,  Patuxent  River,  Maryland 


UNIVERSITIES 


Carleton  University,  Ottawa,  Ont. 

Cornell  University,  Ithaca,  N.Y. 

Ecole  Polytcchnique,  Montreal,  P,Q. 

Laval  University,  Quebec,  P.Q. 

McGill  University,  Montreal,  P.Q. 

Queen's  University,  Kingston,  Ont, 

State  University  of  New  York,  Albany,  N.Y. U.S.A. 
University  of  Alberta,  Edmonton,  Alta. 

University  of  British  Columbia,  Vancouver,  B.C. 
University  of  Ottawa,  Ottawa,  Ont. 

University  of  Sherbrooke  Sherbrooke,  P.Q.: 
University  of  Surrey,  Surrey,  England 
University  of  Toronto,  Toronto,  Ont. 

University  of  Waterloo,  Waterloo,  Ont. 


